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20.  ABSTRACT  (Continued) 


> weapon  burst  height  may  vary  from  surface  to  a maximum  of 
1000  meters  above  ground  and  yields  from  0.01  to  500  kilo- 
tons  were  considered.  The  calculational  methodology  is 
described  and  the  sensitivity  of  the  results  to  various 
quantities  is  presented.  Results  for  selected  source 
spectra  are  given  and  compared  with  results  from  the 
literature.  Conclusions  concerning  the  data  base  are 
presented . 
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1.  INTRODUCTION 


The  primary  or  initial  radiation  from  the  detonation 
of  a nuclear  weapon,  consists  of  prompt  and  delayed  compo- 
nents. The  prompt  radiation  is  composed  of  the  neutrons 
and  gamma  rays  from  the  fission  and/or  fusion  process  and 
the  secondary  gamma  rays  resulting  from  neutron  capture  in 
the  air  and  ground.  This  component  is  a function  of  the 
weapon  source  spectrum  and  the  burst  height  relative  to  the 
ground.  The  delayed  radiation  is  composed  of  the  neutrons 
and  gamma  rays  emitted  by  the  decay  of  fission  products  in 
the  fireball  as  well  as  the  secondary  gamma  rays  produced 
by  capture  of  the  delayed  neutrons  in  the  air  and  ground. 

This  component  of  the  radiation  is  also  a function  of  the 
burst  height.  Because  of  the  perturbed  environment  produced 
by  the  shock  wave,  the  delayed  component  is  also  a function 
of  the  blast  yield  of  the  weapon.  The  response  of  a detector 
to  these  two  radiation  components  depends  on  the  detectc- 
response  as  a function  of  incident  particle  energy  and  the 
detector  position  relative  to  the  burst  and  ground. 

The  purpose  of  this  work  is  to  develop  a data  base 
which  can  be  used  for  predicting  the  dose  and  the  neutron 
to  gamma  ray  dose  ratio  due  to  the  initial  radiation  from 
the  detonation  of  a nuclear  weapon  in  an  air  over  ground 
geometry.  The  data  are  to  be  applicable  for  the  analysis 
of  tactical  nuclear  weapons  and  thus  imply  low  to  inter- 
mediate yield  weapon  bursts  near  the  ground  with  the  dose 
determined  in  a plane  near  the  air-ground  interface.  There- 
fore, the  weapon  yields  considered  range  from  0.01  to  500 
kilotons  with  burst  heights  ranging  from  surface  to  1000 
meters  above  the  ground.  The  detector  plane  is  located 
1.75  meters  above  the  ground  and  extends  out  to  a radius 
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of  1400  meters.  The  two  responses  of  interest  are  free 
field  and  armor  shielded  (representing  personnel  within 
an  armored  vehicle)  tissue  dose  detectors.  The  specific 
dose  quantities  of  interest  for  both  the  prompt  and  delayed 
components  consist  of  the  following: 

1)  Free-field  neutron  tissue  dose; 

2)  Free-field  gamma  ray  tissue  dose; 

3)  Free-field  secondary  gamma  ray  tissue  dose 
from  neutron  capture  in  the  air  and  ground; 

4)  Armor  shielded  neutron  tissue  dose; 

5)  Armor  shielded  gamma  ray  tissue  dose; 

6)  Armor  shielded  secondary  gamma  ray  tissue  dose 
from  neutron  capture  in  the  air  and  ground;  and 

7)  Armor  shielded  secondary  gamma  ray  tissue  dose 
from  neutron  capture  in  the  armor. 

The  techniques  used  in  generating  the  data  base  are 
outlined  in  the  Methodology  section  of  this  report.  A 
Preliminary  Investigation  section  presents  a brief  cross 
section  sensitivity  analysis  together  with  basic  calcu- 
lational  assumptions  and  pertinent  fundamental  data.  The 
Calculations  section  presents  the  models  employed  in  the 
prediction  of  the  prompt  delayed  radiation  components  and 
also  discusses  various  aspects  of  the  calculations.  Several 
forms  of  the  results  for  a few  selected  weapon  burst 
scenarios  are  shown  in  the  Results  section.  A Comparisons 
section  presents  comparisons  of  the  present  work  with 
results  obtained  from  the  literature  and  from  other  calcu- 
lations. The  Summary  section  discusses  some  of  the  conclu- 
sions concerning  the  results  of  this  work. 
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2.  METHODOLOGY 


It  is  generally  recognized  that  to  adequately  predict 
the  transport  of  radiation  in  an  air  over  ground  geometry 
for  the  burst-detector  geometries  of  interest,  the  effect 
of  the  ground  must  be  included.  This  effect  on  the  prompt 
radiation  can  be  quite  important,  especially  for  predicting 
the  dose  due  to  secondary  gamma  rays.  For  a specified 
weapon  blast  yield,  the  ground  affects  the  shock  wave  and 
the  rate  of  fireball  rise,  both  of  which  influence  the 
time  integrated  delayed  dose. 

2.1  PROMPT  RADIATION 

Since  the  dose  from  any  source  spectrum  is  required, 
it  is  convenient  to  solve  the  adjoint  form  of  the  transport 
equation.  The  detector  response  as  a function  of  energy 
is  then  the  source  term  and  the  transport  equation  is  solved 
for  the  adjoint  flux,  also  referred  to  as  the  source  energy 
importance.  The  source  energy  importance  represents  the 
importance  or  value  of  a source  particle  of  a given  energy 
at  a given  source-detector  distance  in  contributing  to  the 
dose  and  forms  the  basis  for  a data  base  to  be  used  in 
subsequent  studies.  Folding  a source  spectrum  with  this 
importance  function  will  give  the  dose.  The  data  base  is 
to  be  used  to  predict  the  dose  for  any  height  of  burst  up 
to  1000  meters  and  for  detector  ground  ranges  from  0.0  to 
1400.0  meters.  Monte  Carlo  methods  for  solving  the  adjoint 
form  of  the  transport  equation  are  available  but  are  judged 
to  be  impractical  in  view  of  the  requirement  that  the  source 
energy  importance  function  is  to  be  defined  at  a large  number 
of  positions.  Two-dimensional  discrete  ordinates  solutions 
to  the  adjoint  transport  equation  satisfy  the  requirements. 
Therefore,  the  two-dimensional  discrete  ordinates  code, 
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DOTSAI ^ ^ , was  used  to  calculate  the  source  energy 
importance  distributions  for  the  two  detector  responses. 

Since  an  analytic  first  collision  source  treatment 
was  utilized,  the  results  of  the  calculations  are  eleven 
matrices  consisting  of  the  importance  of  uncollided  and 
scattered  source  particles  in  contributing  to  the  dose 
type  as  a function  of  source  particle  energy  and  spatial 
distribution.  These  eleven  distributions  are  for  the 
importance  of: 

1)  Uncollided  neutrons  in  contributing  to  the 
neutron  free-field  tissue  dose; 

2)  Scattered  neutrons  in  contributing  to  the 
neutron  free-field  tissue  dose; 

3)  Scattered  neutrons  in  contributing  to  the  air 
and  ground  induced  secondary  gamma  ray  free- 
field  tissue  dose; 

A)  Uncollided  gamma  rays  in  contributing  to  the 
gamma  ray  free-field  tissue  dose; 

5)  Scattered  gamma  rays  in  contributing  to  the 
gamma  ray  free-field  tissue  dose; 

6)  Uncollided  neutrons  in  contributing  to  the 
neutron  armor  shielded  tissue  dose; 

7)  Scattered  neutrons  in  contributing  to  the 
neutron  armor  shielded  tissue  dose; 

8)  Scattered  neutrons  in  contributing  to  the  air 
and  ground  induced  secondary  gamma  ray  armor 
shielded  tissue  dose; 
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9)  Scattered  neutrons  in  contributing  to  the  armor 
induced  secondary  gamma  ray  armor  shielded 
tissue  dose; 

10)  Uncollided  gamma  rays  in  contributing  to  the 
gamma  ray  armor  shielded  tissue  dose;  and 

11)  Scattered  gamma  rays  in  contributing  to  the 
gamma  ray  armor  shielded  tissue  dose. 

2.2  DELAYED  RADIATION 

The  solution  technique  employed  for  the  prompt 
radiation  component  is  impractical  for  the  delayed  rad- 
iation. The  fission  products  are  in  the  rising  and  ex- 
panding fireball  and  thus  the  geometry  changes  with  time. 
The  time  decay  of  the  source  spectra  for  the  delayed 
components  is  assumed  to  be  independent  of  energy.  The 
air  density  is  determined  by  the  shock  wave  and  the  source 
intensity  viewed  by  the  detector  is  affected  by  the  rate 
of  fireball  rise.  Both  of  these  properties  are  functions 
of  the  weapon  blast  yield  and  burst  height.  Thus,  for  a 

given  burst  scenario  the  delayed  component  scales  directly 
only  with  the  fission  yield.  Because  of  these  factors  a 
different  technical  approach  was  utilized  for  the  delayed 
radiation . 

Existing  methods  for  predicting  dose  from  the  delayed 
radiation  component  are  based  on  the  assumption  that  4ttr2 
dose  in  air  scales  directly  with  rho-R , often  referred  to 
as  rho-R  scaling,  where  rho-R  is  the  mass  thickness  of  air 
between  the  source  and  the  detector.  Obviously,  any 
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calculational  method  must  model  the  shock  wave  - fireball 

( o 3) 

phenomenology.  The  IDEA'  ’ code  incorporates  such  a 
model  and  is  designed  to  predict  the  prompt  and  delayed 
integrated  dose  in  air.  The  shock  wave  and  fireball 
phenomenology  include  the  effect  of  the  ground.  Therefore, 
the  IDEA  code  was  utilized  in  the  following  manner.  The 
delayed  neutron  and  fission  product  gamma  ray  source  spectra 
for  a unit  source  strength  were  folded  with  the  source 
energy  importance  distributions,  resulting  in  the  delayed 
dose/source  particle  in  uniform  density  air  for  all  burst 
heights  and  detector  locations  of  interest.  The  various 
components  were  then  converted  to  4nR  dose  as  a function 
of  rho-R.  As  the  time  incremented  blast  phenomenology  is 
treated  in  IDEA,  dose  rates  are  derived  from  this  matrix 
utilizing  rho-R  scaling.  The  advantages  of  this  procedure 
are  that  the  effect  of  the  ground  on  the  dose  is  thus 
implicit  and  a consistent  cross  section  data  base  is  utilized. 
The  disadvantage  is  that  many  weapon  burst  scenarios  must 
be  investigated  to  generate  the  appropriate  delayed  dose 
portion  of  the  data  base.  Techniques  must  then  be  developed 
for  applying  the  results  to  environments  not  explicitly 
calculated. 

In  summary,  the  fundamental  calculational  technique 
is  the  solution  of  the  adjoint  form  of  the  transport  equation 
to  determine  source  particle  importance  distributions  as  a 
function  of  source  energy  and  source  location.  The  imple- 
mentation of  this  method  for  predicting  both  the  prompt 
and  delayed  dose  components  will  be  described  in  more 
detail  in  the  Calculations  section. 
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3.  PRELIMINARY  INVESTIGATIONS 


The  first  phase  of  the  effort  to  prepare  a tactical 
nuclear  weapons  radiation  environment  data  base  consisted 
of  preliminary  investigations  to  determine  an  appropriate 
cross  section  data  base  and  group  structure,  the  sensitivity 
of  results  to  ground  composition,  appropriate  response 
functions,  etc.  Each  of  these  items  as  well  as  other  input 
data  are  discussed  in  this  section. 

3.1  CROSS  SECTIONS 

Previous  calculations  of  radiation  environments  for 
an  air  over  ground  geometry  were  performed  using  cross 
sections  of  the  1967  era.  More  recent  point  cross  section 
data  have  been  incorporated  in  the  Evaulated  Nuclear  Data 
File,  Version  IV  (ENDF/IV)  used  in  processing  a DNA  coupled 
37-21  neutron-gamma  ray  library  maintained  and  distributed 
by  the  Radiation  Shielding  Information  Center  (RSIC)  as 
DPL-31.  Another  multigroup  library  generated  in  1970  has 
been  distributed  by  RSIC  as  the  CASK  Library  (DLC-23) . 
Additionally,  a 22-18  neutron-gamma  ray  SAI  library  was 
processed  in  the  Spring  of  1974  from  the  best  point  data 
available  including  the  ENDF/IV  data. 

3.1.1  Air  Cross  Sections 

The  transport  of  radiation  in  an  air  over  ground 
geometry  is  dominated  by  the  cross  sections  of  nitrogen 
and  oxygen  with  the  importance  of  the  ground  cross  sections 
appearing  in  the  calculation  of  the  secondary  gamma  ray 
component.  To  determine  the  adequacy  of  the  nitrogen  cross 
section  - the  dominant  element  for  air  - comparisions  of 
calculations  with  experiments  were  made.  Previous  work ^ has 
demonstrated  the  differences  to  be  expected  in  results  calcu- 
lated with  the  older  cross  sections  (frequently  referred  to  as 
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the  Straker  set) and  the  newer  ones  (frequently  referred 
to  as  the  Young  set)  . 

The  agreement  between  experiment  and  calculation  is 
used  as  a measure  of  the  state-of-knowledge  of  the  cross 
sections.  Recent  integral  experiments  at  Intelcom  Rad- 
iation Tech  (IRT)^)  and  Christian  Albrechts  University^) 
offer  new  results  which  can  be  used  to  evaluate  the  adequacy 
of  differential  nuclear  data. 

3 . 1.1.1  1RT  Pulsed  Sphere  Experiment 

(9) 

The  MORSE  code  with  combinatorial  geometry  was  used 
to  calculate  the  IRT  pulsed  sphere  experiment.  The  exper- 
imental geometry  was  modeled  by  including  only  the  spherical 
styrofoam  container  and  the  liquid  nitrogen.  Both  forward 
and  adjoint  calculations  were  performed  using  the  Straker 
and  Young  cross  section  sets.  The  experimental  results 
are  presented  as  a function  of  incident  neutron  energy  and 
therefore  adjoint  transport  calculations  can  be  directly 
compared  with  the  experimental  results  and  are  independent 
of  the  neutron  source.  The  forward  calculations  were  used 
to  more  efficiently  investigate  the  gamma  ray  counts  due 
to  neutrons  incident  on  the  detector. 

The  detector  efficiencies  given  in  Table  1 were  used 
as  response  functions  in  the  forward  calculations  and  as 
sources  in  the  adjoint  calculations.  The  gamma  ray  response 
for  neutrons  incident  on  the  detector  is  included  in  the 
detector  efficiency  for  gamma  rays. 

In  the  forward  calculations  the  detector  counts  were 
scored  as  a function  of  time.  The  time  and  energy  distri- 
bution of  source  neutrons  was  utilized.  In  order  to  compare 
with  the  adjoint  results,  the  detector  counts  were  trans- 
formed from  time  space  to  energy  space  and  divided  by  the 

source  intensity  to  give  detector  counts  per  incident 

2 

neutron/cm  versus  incident  neutron  energy. 


i 
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Table  1.  Detector  efficiency. 


Inc ident 

Energy 

Particle 

(MeV) 

Neutrons 

12.2 

- 

15.0 

10.0 

- 

12.2 

8. 18 

- 

10.0 

6.36 

- 

8. 18 

4.96 

- 

6.36 

4.06 

- 

4.96 

3 01 

- 

4.06 

2.46 

- 

3.01 

2.35 

- 

2.46 

1.826 

- 

2.35 

1.108 

- 

1.  826 

0.55 

- 

1.108 

Gamma  Rays 

8.0 

- 

10.0 

6.5 

- 

8.0 

5.0 

- 

6.5 

4.0 

- 

5.0 

3.0 

- 

4.0 

2.  5 

- 

3.0 

2.0 

- 

2.5 

1.66 

- 

2.0 

1.33 

- 

1.66 

1.0 

- 

1.33 

0.8 

- 

1.0 

0.6 

- 

0.8 

0.4 

- 

0.6 

0.3 

- 

0.4 

Counts  Per  Incident  Particle 

Neutron 

Gamma  Ray 

0.  142 

0.  0116 

0.  150 

0. 0134 

0.  162 

0.0125 

0.  180 

0.  01 

0.208 

0.  0075 

0.203 

0.  00525 

0.249 

0.  00375 

0.  266 

0.  0025 

0.270 

0.  0019 

0.270 

0.  0016 

0.244 

0.  00125 

0.  060 

0.  001 

0.0 

0.  16 

0.0 

0.  163 

0.0 

0.  166 

0.0 

0.  170 

0.0 

0.  173 

0.0 

0.  176 

0.0 

0.  179 

0.0 

0.  183 

0.0 

0.  185 

0.0 

0.  186 

0.0 

0.  182 

0.0 

0.  170 

0.0 

0. 150 

0.0 

0.  070 
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The  results  of  the  adjoint  transport  calculations 
for  the  detector  efficiencies  are  presented  for  neutron  and 
gamma  ray  counts.  The  calculations  are  compared  with  the 
experimental  results  for  a detector  at  30°  in  Figure  1.  For 
the  neutron  results  at  the  30°  detector  position,  the  Young 
cross  section  set  gives  results  which  are  15  - 20  percent  lower 
than  the  measurements  in  the  energy  range  below  7 MeV.  The 
Straker  cross  sections  given  even  poorer  results.  The 
results  of  the  secondary  gamma  ray  calculations  using  the 
Young  cross  sections  give  good  agreement  (within  10  percent) 
with  the  measurements  for  all  of  the  detector  positions  and 
for  neutron  energies  greater  than  2-3  MeV.  The  Straker 
cross  sections  give  results  which  typically  are  a factor 
of  2 high  for  incident  neutron  energies  above  9 MeV  and  lower 
than  the  measurements  for  energies  below  8 MeV. 

3. 1.1.2  Christian  Albrechts  University  Tank  Experiment 

This  tank  experiment  is  similar  in  many  respects  to 

the  GRT  liquid  nitrogen  dewar  exper iment which  has 

previously  been  calculated  using  the  Straker  and  Young  cross 

section  sets  for  nitrogen  . Similar  calculations  have 

been  performed  for  the  Christian  Albrechts  University  (CAU) 

tank  experiment.  One  dimensional  multigroup  discrete 

ordinates  calculations  were  performed  with  the  ANISN  rad- 

112) 

iation  transport  code  . The  calculations  were  made  with 

a unit  source  in  the  highest  energy  group.  For  comparison 

with  the  published  experimental  results,  the  neutron  energy 

2 

fluxes  were  converted  to  4tiR  fluence  per  unit  lethargy  per 

source  neutron.  The  secondary  gamma  ray  intensities  are 
2 

reported  as  4ttR  photons  per  MeV  per  source  neutron. 

A comparison  of  the  results  of  the  tank  experiment  with 
calculations  using  Young  and  Straker  cross  sections  for  neu- 
trons and  secondary  gamma  rays  is  shown  in  Figures  2 and  3, 
respectively.  For  neutrons,  there  is  better  agreement 
between  measurements  and  calculations  using  Young  cross 
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Figure  1.  Calculated  and  measured  scattered  neutron  and 


secondary  gamma  ray  spectra  at  30°  from 
Nitrogen  Sphere. 
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Figure  2.  Calculated  and  measured  neutron  spectrum  at 
60  cm  from  a 14  MeV  neutron  source  in 
liquid  nitrogen. 
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Garrma  Ray  Intensity  (Photon*/ McV/Source  Neutron) 


sections  than  when  Straker  cross  sections  are  used.  Regions 
of  disagreement  are  in  the  source  energy  region  around  15  MeV 
and  around  1.5  MeV  where  the  experiment  is  higher.  There  are 
fairly  minor  differences  between  the  two  calculations  except 
in  the  region  of  the  source  energy. 

For  secondary  gamma  rays,  there  are  significant  differ- 
erences  between  the  calculations  and  measurements.  Both 
calculations  fall  substantially  below  the  experimental  meas- 
urements with  the  calculations  based  on  Young  cross  sections 
being  about  a factor  of  2 lower  than  the  measurement.  Although 
the  agreement  is  better  when  Straker  cross  sections  are  used, 
the  shape  of  the  calculated  curve  does  not  agree  with  the 
measurement.  One  possibility  of  the  disagreement  is  that 
gamma  ray  production  in  the  detector  is  significant  at  the 
deeper  penetration  distances  and  this  effect  was  not  considered 
in  the  reduction  of  the  experimental  data  and  could  not  be 
included  in  the  analysis  since  the  detector  response  was  not 
provided. 

In  summary,  calculations  based  on  Young  cross  sections 
predict  the  IRT  experiment,  whereas  calculations  based  on 
Straker  cross  sections  provide  results  which  agree  better 
with  the  CAU  tank  experiment.  Due  to  the  improved  data  on 
which  the  Young  evaluation  was  based,  the  Young  evaluated  set 
was  used  for  the  data  base  calculations. 

3.1.2  Ground  Cross  Sections  and  Composition 

There  are  two  components  that  affect  the  ground 
contribution  to  the  radiation  environment.  These  are  the 
ground  composition  and  the  ground  cross  sections  and  both 
of  these  factors  were  investigated. 

A typical  soil  is  desired  in  the  sense  that  its 
composition  affects  the  radiation  environment.  However, 
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the  ground  composition  to  he  used  in  the  calculations  should 
te  selected  with  some  justification  in  order  for  the  results 
t"  he  generally  useful.  An  analysis  of  the  effects  of  ground 
composition  and  moisture  content  on  the  radiation  environments 
produced  at  the  air-ground  interface  was  performed.  The 
calculat ional  model  consisted  of  a small  yield  thermonuclear 
weapon  denotated  at  a point  100  meters  above  the  ground 
surface.  Consideration  was  given  to  both  a tissue  dose  and 
an  armor  shielded  tissue  dose  response  at  a point  50  centi- 
meters above  the  ground  surface. 

At  the  time  of  the  analysis,  sensitivity  calculations 
were  restricted  to  one  dimensional  geometry.  The  best  one 
dimensional  approximation  of  the  air  over  ground  problem 
for  the  investigation  of  ground  effects  is  a slab  geometry 
representation.  Figure  4 shows  the  calculational  model 
used  in  the  one  dimensional  transport  calculations.  An  air 
albedo  boundary  condition  was  used  above  700  meters.  Only 
the  first  50  centimeters  of  ground  was  included  in  the 
calculation  and  a vacuum  boundary  condition  was  used  at  the 
bottom  of  the  problem.  Both  of  these  boundary  conditions 
are  felt  to  be  adequate. 

The  calculational  approach  was  to  assume  a baseline 
soil  composition  to  investigate  the  effects  of  varying  the 
composition  using  first  order  perturbation  theory.  First 
order  perturbation  theory  assumes  that  the  effects  of  changes 
in  soil  composition  on  the  free -field  and  armor  shielded 
tissue  dose  are  small  (less  than  10-15  percent).  The  sensi- 
tivity of  the  air-ground  radiation  environment  was  determined 
individually  for  a number  of  elements  and  the  predicted 
perturbations  due  to  various  soils  determined  from  the  element 
sensitivities.  Forward  and  adjoint  calculations  were  performed 
with  the  one  dimensional  discrete  ordinates  transport  code  ANISN. 
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4.  One  dimensional  calculational  model  used 
for  sensitivity  analysis. 
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A first  order  perturbation  analysis  was  performed  using  the 

ANISN  results  and  the  one  dimensional  cross  section  sensi- 

(in 

tivity  analysis  code,  SWANLAKE  . One  of  the  results 
obtained  from  SWANLAKE  is  the  "sensitivity"  due  to  the 
perturbation  cross  section  set.  The  "sensitivity"  is  equal 
to  the  fractional  change  in  the  total  result  due  to  the 
perturbation  cross  section,  AE.  That  is,  the  sensitivity, 


(1) 


where 

R is  the  total  response,  and 

AR  is  the  perturbation  to  the  total  response. 


If  the  sensitivities,  S^,  are  determined  for  a number  of 

elements,  i,  with  densitv  1.- , the  sensitivity  of  a 

7 barn* cm  ’ 7 

mixture  of  elements,  S . , is  given  by 

mi  V ° J 


# elements 


mix 


£ 


N.S. 


i i . 


(2) 


where 


are  the  element  densities  in  the  mixture. 

In  order  to  estimate  the  perturbation  to  the  total  response 
due  to  a different  soil  composition,  equations  (1)  and  (2) 
can  be  combined  to  give, 


# elements 
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where 

R is  the  perturbed  response 

is  the  element  density  for  the  baseline  soil,  and 
is  the  element  density  for  the  desired  soil. 

The  results  of  the  forward  and  adjoint  ANISN  calcula- 
tions for  the  free-field  and  armor  shielded  tissue  dose  are 
given  in  Table  2.  The  agreement  of  the  forward  and  adjoint 
calculations  is  better  than  0.57,.  The  armor  shielded  res- 
ponse due  to  incident  gamma  rays  is  reduced  by  approximately 
a factor  of  4 for  the  shielded  case.  For  this  calculational 
model  approximately  907,  of  the  free-field  response  is  due  to 
incident  neutrons  while  977«  of  the  armor  shielded  response 
is  due  to  incident  neutrons.  The  ratio  of  dose  components 
is  strongly  dependent  on  the  source-detector  range  and  thus 
is  different  for  larger  distances. 

Several  soil  compositions  were  considered.  Table  3 gives 
the  atomic  concentrations  (^>15)  £or  a number  Gf  sofis  found 
in  the  United  States  and  Western  Europe.  Results  for  the 
sensitivity  of  ground  composition  for  several  elements  are 

given  in  Table  4 for  both  the  free-field  and  armor  shielded 
response.  It  is  noted  that  hydrogen  is  the  most  important 
element  from  the  sensitivty  analysis.  Other  elements  such 
as  titanium,  manganese  and  molybdenum  are  also  important. 

The  hydrogen  content  can  be  quite  large  in  soils  and  vary 
substantially  due  to  the  moisture  content.  The  hydrogen 
sensitivity  is  negative  which  means  that  perturbation  theory 
predicts  that  the  tissue  dose  would  decrease  for  a small 
increase  in  hydrogen  density.  For  the  free-field  case  most 
of  the  large  positive  sensitivities  result  from  secondary 
gamma  ray  production  by  these  elements  in  the  soil.  However, 
for  the  shielded  tissue  response  t ne  secondary  gamma  rays 
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Table  2.  Forward  and  adjoint  calculated 
baseline  tissue  dose. 


Free-Field 

Total 

Neutron 
Gamma  Ray 
Armor  Shielded 
Total 

Neutron 
Gamma  Ray 


Tissue  Dose  50  cm  Above  Ground 
(rads /source  neutron) 

Forward 

Adjoint 

2.  682  x 10"9 

2. 688  x 10"9 

2.412  x 10'9 

2.702  x 10"10 

2.496  x 10'9 

2.500  x 10'9 

2.430  x 10"9 

6. 600  x 10'11 

Table  3.  Composition  of  United  States  and  Western  European  Soils. 


Table  4.  Sensitivity  of  calculated  total  dose 
to  ground  composition. 


Element 

Sensitivity3 

Element 

2^ 

Sensitivity 

Free-Field 

Armor 

Shielded 

Free-Field 

Armor 

Shielded 

H 

- 33.  0 

- 26.5 

K 

5.24 

- 8.40 

C 

0.  595 

1.65 

Ca 

0.0149 

- 1.04 

0 

l 34 

2.24 

Ti 

26.6 

- 13.  6 

Na 

0.462 

0.  892 

Cr 

0.  798 

- 7.32 

Mg 

1.93 

2.  52 

Mn 

18.0 

- 50.2 

A1 

2.25 

1.48 

Fe 

8.  63 

- 5.4 

Si 

1.33 

1.34 

Mo 

13.5 

- 1. 66 

‘ Fractional  change  in  dose  due  to  a 1.0  atom  (barn  cm)  addition  of  the 
element. 
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from  the  ground  are  not  as  important  and  the  sensitivity  is 
in  fact  often  negative.  Although  the  sensitivity  calculations 
were  performed  for  a simple  one-dimensional  model,  the  results 
indicate  that  the  tissue  dose  may  vary  by  a factor  of  two  or 
more  depending  on  the  hydrogen  content  of  the  soil. 

The  effect  of  moisture  in  the  soil  can  be  estimated 
from  the  element  sensitivities  for  hydrogen  and  oxygen.  In 
this  manner  a unit  weight  percent  increase  in  moisture  content 
would  result  in  a -3.47,  and  -2.77,  change  in  the  free-field  and 
armor  shielded  tissue  dose,  respectively.  This  predicted 
moisture  content  sensitivity  was  checked  by  a direct  recalcu- 
lation of  the  baseline  problem  with  the  addition  of  1 weight 
percent  water  and  the  agreement  with  the  perturbation  calcu- 
lation was  quite  good. 

The  element  sensitivities  were  used  to  estimate  the 
correction  or  perturbation  factor  given  in  equation  3 for 
all  of  the  soil  types  listed  in  Table  2.  The  results  are 
given  in  Table  5.  In  Figure  5,  the  ratio  of  the  dose  is 
plotted  as  a function  of  the  hydrogen  content  of  the  various 
soils.  It  is  observed  that  most  of  the  variations  in  the 
dose  ratio  can  be  explained  by  the  hydrogen  content.  The 
maximum  differences  observed  ranged  from  +157,  to  -387,. 

From  this  investigation  which  utilized  a one-dimensional 

model,  it  was  concluded  that  an  area  weighted  average  Western 

European  soil  composition  was  a reasonable  ground  composition 

to  use  in  performing  the  calculations.  The  air-ground  tissue 

response  is  fairly  sensitive  to  the  hydrogen  content  of  the 

_ 3 

soil  and  a hydrogen  density  greater  than  2.0  x 10  atoms/ 
(barn  • cm)  results  in  a significantly  reduced  dose.  The 
compositions  of  ground  and  air  used  in  the  data  base  calcu- 
lations are  given  in  Table  6. 
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Table  5.  Sensitivity  of  calculated  total  dose 
to  soil  type. 


Correction  Factor,  Ca 

Soil  Type 

Free  Field 

Armor  Shielded 

Mean  Earth's  Crust 

1.0 

1.0 

Central  German  Soils 

sand 

1.03 

1.06 

loess 

0.99 

1.02 

clay 

0.  71 

0.  78 

marl  loam 

0.92 

0.95 

shale 

0.  84 

0.  88 

sandstone 

0.98 

1.0 

marl 

0.  84 

0.  88 

limestone 

1.01 

1.04 

pranite 

1.  03 

1. 04 

bassalt 

0.99 

0.98 

pneiss 

0.61 

0.  70 

slate 

0.92 

0.95 

tsoil 

1.  15 

0.  85 

area  weiphted 

0.98 

1.  01 

United  States  Soils 

beach  sand 

1.05 

1.  08 

lava  clay 

0.72 

0.72 

Nevada  dessert 

0.  82 

0.  82 

The  correction  factor  is  the  ratio  of  the  dose  expected  for  the 
particular  soil  and  the  dose  obtained  for  the  baseline  case. 


% 
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RATIO  OF  DOSE  TO  BASE  CASE  VALUE 


Table  6.  Air  and  ground  composition  for  tactical  nuclear 
weapons  data  base  calculations. 


Element 

Atoms  (barn  • centimeter) 

A a 

Air 

(1.11  x 10"^  gm /cm**) 

Ground 

3 

(1.6  gm  c m ) 

H 

1.  753  (-3) 

C 

1. 639  (-3) 

N 

3.635  (-5) 

O 

9.620  (-6) 

3.035  (-2) 

Na 

1.586  (-4) 

Mr 

2.577  (-4) 

A1 

6.969  (-3) 

Si 

1.  093  (-2) 

K 

4.922  (-4) 

Ca 

1.529  (-3) 

Fe 

2.795  (-4) 

aFinal  data  base  results  rho-R  scaled  to  1.  19  mg  cm'*  air  density  (see 
Section  4.2) 


The  other  component'  of  the  effect  of  the  ground  is  that 
due  to  the  cross  sections  for  the  elements  in  the  ground.  The 
same  calculational  model  used  to  investigate  the  effect  of 
composition  was  also  used  to  determine  the  effect  of  using 
different  cross  section  sets.  The  CASK,  SAI  and  DNA  libraries 
were  used  to  calculate  the  free  field  and  armor  shielded 
tissue  dose  due  to  a thermonuclear  source.  Due  to  the  lack  of 
cross  sections  for  some  elements  in  the  libraries,  different 
ground  compositions  were  used  and  are  given  in  Table  7.  The 
results  of  the  calculations  are  shown  in  Table  8.  It  is  noted 
that  the  different  cross  sections  have  little  effect  (less 
than  570)  on  the  results. 

The  cross  section  set  used  for  the  calculations  was  the 
DNA  library.  In  order  to  reduce  computational  effort  and  the 
size  of  the  resulting  data  base,  the  cross  sections  were  col- 
lapsed with  ANISN  to  a 23-19  neutron-gamma  ray  library  using 
the  spectrum  in  air  due  to  a thermonuclear  source.  The  group 
structure  is  given  in  Table  9. 


Table  7.  Soil  composition  used  with  each  cross 
section  sec. 


Element 

Soil  Composition  (atoms/barn  • cm) 

CASK  22-18 
Average  Soil 

SAI  22-18 
Average  Soil 

DNA  38-21 
Average  Soil 

H 

1.753  (-3) 

1.753  (-3) 

1.753  (-3) 

C 

1.  639  (-3) 

1.  639  (-3) 

1.  639  (-3) 

O 

3.  035  (-2) 

3.035  (-2) 

3.035  (-2) 

Na 

1.586  (-4) 

0.0 

1.  586  (-4) 

Mr 

2.577  (-4) 

2.577  (-4) 

2.577  (-4) 

A1 

6.969  (-3) 

6.969  (-3) 

6.969  (-3) 

Si 

1.  093  (-2) 

1.093  (-2) 

1.  093  (-2) 

K 

4.922  (-4) 

0.  0 

4.922  (-4) 

Ca 

1.529  (-3) 

1.529  (-3) 

1.529  (-3) 

Ti 

1.411  (-5) 

1.411  (-5) 

0.  0 

Mn 

1. 833  (-6) 

0.0 

0.0 

Fe 

2.795  (-4) 

2.795  (-4) 

2.  795  (-4) 

Table  8.  Calculated  total  tissue  dose  as  a 
function  of  cross  section  set. 


Cross  Section  Set 


CASK  22-18 


SAI  22-18 


DNA  38-21 


Total  Tissue  Dose  (rads  source  neutron) 


Shielded 


2.  529  (-9) 


2.493  (-9) 


2.  600  (-9) 


Unshielded 


2.  661  (-9) 


2.593  (-9) 


2.  715  (-9) 


Table  9. 


Energy  group  structure  ol  tactical  nuclear 
weapons  radiation  environment  data  base. 


Neutron 

Group 

— 

Energy 

(eV) 

Gamma  Ray 
Group 

Energy 

(eV) 

1 

1.  5 (+7)  - 1.22  (+7) 

1 

1.0  ( + 7)  - 8.0  (+6) 

2 

1.22  (+7)  - 1.00  (+7) 

2 

8.0  (+6)  - 6.0  (+6) 

3 

1.00  (+7)  - 8.18  (+6) 

3 

6.0  (+6)  - 5.0  (+6) 

4 

8.  18  (+6)  - 6.37  (+6) 

4 

5.0  (+6)  - 4.0  (+6) 

5 

6.37  (+6)  - 4.96  (+6) 

5 

4.0  (+6)  - 3.0  (+6) 

6 

4.96  (+6)  - 4.06  (+6) 

6 

3.0  (+6)  - 2.5  (46) 

7 

4.06  (+6)  - 3.01  (+6) 

7 

2.5  (+6)  - 2.0  (+6) 

8 

3.01  (+6)  - 2.38  (+6) 

8 

2.0  (+6)  - 1.5  (+6) 

9 

2.38  (+6)  - 2.30  (+6) 

9 

1.5  (+6)  - 1.0  (+6)  1 

10 

2.30  (+61  - 1.83  (+6) 

10 

1.0  (+6)  - 7.0  (+5) 

11 

1.83  (+6)  - 1.11  (46) 

11 

7.0  (+5)  - 4.5  (+5) 

12 

1.11  (+6)  - 5.50  (+5) 

12 

4.5  (+5)  - 3.0  (+5)  J 

13 

5.50  (+5)  - 1.  11  (+5) 

13 

3.0  (+5)  - 1.5  (+5) 

14 

1.  11  (^5)  - 2.  18  ( + 4) 

14 

1.5  (+5)  - 1.0  (+5)  , 

15 

2. 18  (.4)  - 3.35  (+3) 

15 

1.0  (+5)  - 7.0  (+4)  | 

16 

3.35  (43)  - 5.83  (+2) 

16 

7.0  (+4)  - 4.5  (+4) 

17 

5.  83  (+2)  - 1.01  (42) 

17 

4.5  (+4)  - 3.0  (+4) 

18 

1.01  (42)  - 2.90  (4I) 

18 

3.0  (+4)  - 2.0  (44) 

19 

2.90  (-1)  - 1.07  (+1) 

19 

20 

1.07  ( + 1)  - 3.06  (4O) 

21 

3.  06  (+0)  - 1.12  (+0) 

22 

1.  12  (4O)  - 4.  14  (-1) 

23 

4.  14  (-1)  - 1.00  (-5) 
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3.2  RESPONSE  FUNCTIONS 


The  biological  response  to  radiation  is  the  quantity  of 
interest  in  this  study.  Tissue  kerma  for  the  free  field,  i.e., 
radiation  environment  in  a flat  earth  geometry,  was  chosen  as 
the  response  function.  It  was  also  desirable  to  provide 
information  on  the  extent  of  shielding  provided  by  a light 
armored  tank.  The  response  function  of  a complex  vehicle 
depends  upon  the  energy-angular  description  of  the  radiation 
incident  upon  the  tank  as  well  as  the  specific  location  of 
the  detector  within  the  tank.  Because  of  random  orientation 
of  a vehicle  with  respect  to  a burst  and  of  the  several  loca- 
tions within  a tank,  a simplified  model  was  used  to  determine 
an  armor  shielded  tissue  response  function. 

A one  dimensional  calculationai  model,  illustrated  in 
Figure  6,  was  used  with  a tissue  kerma  response  at  the  center. 

The  armor  composition  is  given  in  Table  10.  Adjoint  calculations 
were  performed  to  determine  the  neutron  and  gamma  ray  response 
of  the  shielded  detector.  The  neutron  response  consists  of  the 
neutron  dose  inside  due  to  incident  neutrons  and  a component 
of  the  gamma  ray  dose  inside  due  to  gamma  rays  produced  in  the 
armor  by  the  incident  neutrons.  The  incident  gamma  rays  also 
contribute  to  the  gamma  ray  dose  inside.  Table  11  gives  the 
free-field  and  armor  shielded  detector  tissue  response  functions. 
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ARMOR 


Figure  6.  Armor  shielded  tissue  dose  response  function 
calculational  model 


Table  10.  bight  tank  armor  composition. 


Element 

Weight  Percent 

Atomic  Concentration 
atoms /(barn  • cm) 

Fe 

95.63 

8.097  x 10"2 

C 

0.27 

1.063  x 10"3 

Mn 

1.03 

8. 866  x 10'4 

Si 

0.48 

8.081  x 10"4 

Cr 

1.07 

9.  729  x 10"4 

Ni 

0.97 

7.813  x 10'4 

Mo 

0.55 

2.711  x 10"4 

Table  11.  Free-field  and  armor  shielded 

detector  tissue  response  functions. 


Tissue  Kerma  Factor 

2 

(rads  • cm  incident  particle) 

Group 
( Neutrons] 

Armor 

Shielded 

Neutron 

Gamma 

Rav 

Total 

1 

1. 50  (.7) 

6.36  (-  9) 

4.  30  (-  9) 

3.67  (-10) 

4. 67  (-  9) 

2 

1.22  (*7) 

5.74  (-  9) 

3.  77  (-  9) 

3.96  (-10) 

4.  17  (-  9) 

3 

1.00  (.7) 

5.  1"  (-  9) 

3.51  (-  9) 

4.03  (-10) 

3.91  (-  9) 

4 

8.  18  (*6) 

4.  87  (-  9) 

3.  20  (-  9) 

3.49  (-10) 

3. 55  (-  9) 

5 

6.37  (-6) 

4.51  (-  9) 

3.  02  (-  9) 

2.  76  (-10) 

3.30  (-  9) 

6 

4.96  (♦6) 

4.21  (-  9) 

2 74  (-  9) 

2.  08  (-10) 

2.  94  (-  9) 

7 

4.  06  (*6) 

3.98  (-  9) 

2.  64  (-  9) 

1.27  (-10) 

2.77  (-  9) 

8 

3.01  (*6) 

3.39  9) 

2.47  (-  9) 

5.  74  (-11) 

2.53  (-  9) 

9 

2.38  (+6) 

3.  07  (-  9) 

2.31  (-  9) 

4. 24  (-11) 

2.35  (-  9) 

10 

2.30  (-6) 

3.05  (-  9) 

2.31  (-  9) 

4.07  (-11  1 

2.35  (-  9) 

11 

1 . 83  i 

2.  63  (-  9) 

2.  08  (-  9) 

2.71  (-11) 

2. 11  (-  9) 

12 

1. 11  (-6) 

2.  AC-  1 1) 

1.  84  (-  9) 

1.08  (-11) 

1.  85  (-  9) 

13 

5.50 

1 . 2 : | - l) 

1.21  (-  9) 

7.  87  (-12) 

1.  22  (-  9) 

14 

1.11 

4 00  (-10) 

3.25  (-10) 

2.75  (-111 

3.53  (-10) 

15 

2.  13  t • 4 1 

4.00  (-10) 

3. 25  (-10) 

2.75  (-11) 

3.  53  (-10) 

16 

3.35  (-3i 

1.96  (-11) 

8.  91  (-12) 

8.  19  (-11) 

9.03  (-11)  | 

17 

5.83  < 2 

1.67  (-12) 

1.40  (-12) 

4.25  (-11) 

4 39  (-11)  j 

18 

1.01  (-2) 

. 17  (-12) 

6.27  (-13) 

5.57  (-Hi 

5.  64  6-11) 

19 

2.90  (.1) 

1.  11  (-121 

8.  12  (-13) 

1.04  (-10) 

1.  04  (-10) 

20 

1.07  (-1  i 

1.  62  (-12) 

8.  27  (-13' 

1.37  (-10) 

1.38  (-10) 

21 

3.  06  (.0) 

2.65  (-12) 

7.68  (-13) 

1.95  (-10' 

1.96  (-10) 

22 

1.  12  (-0) 

4.26  (-12) 

5.67  (-13) 

3. 07  (-10) 

3.  08  (-10) 

23 

[ Gamma 

Rays] 

4. 14  (-1) 

9.35  (-12) 

1.  13  (-13) 

6.47  (-10) 

6.47  (-10) 

24 

1. 00  (-7) 

2.42  (-  9) 

0.0 

7.  12  (-10) 

7.  12  6- 1 0) 

25 

8.00  (.6) 

1.95  (-  9) 

0.  0 

6.  03  (-10. 

6.  03  < -10)  , 

26 

6.  00  (.6) 

1. 84  (-  9) 

0.  0 

5.77  (-10) 

5.77  (-10) 

27 

5.00  (*6) 

1.59  (-  9) 

0.0 

4 74  (-10) 

4.  74  (-10) 

28 

4.  00  ( ♦ 6 ) 

1.27  (-  9) 

0.0 

3. 67  (-10) 

3.  67  (-10) 

23 

3.00  (.6) 

1.03  (-  9) 

0.0 

2.  86  (-10) 

2.  86  (-10)  I 

30 

2.50  (46) 

8.75  (-10) 

0.  0 

2.  16  (-10) 

2.  16  (-10) 

31 

2.  00  (+6) 

7.05  (-10) 

0.  0 

1.49  (-10) 

1.49  (-10)  | 

32 

1.50  (<6) 

5.70  (-10) 

0.  0 

9.35  (-11) 

9.35  (-11) 

33 

1.00  ( + 6) 

4. 13  (-10) 

0.0 

4. 18  (-11) 

4.  18  (-11) 

34 

7.  00  (.5) 

2.94  (-10) 

0.  0 

1.91  (-11) 

1.91  (-11) 

35 

4.50  (*5) 

2.  03  (-10) 

0.  0 

7.01  (-12) 

7.01  (-12) 

36 

3.  00  ( ♦ 5 ) 

1.03  (-10) 

0.  0 

7.04  (-131 

7.04  (-13) 

37 

1.50  (*5) 

6.60  (-11) 

0.0 

1.22  (-14) 

1.22  (-14) 

38 

1.00  (.5) 

3.90  (-11) 

0.0 

0.0 

0.  0 

39 

7.  00  (44) 

4.79  (-11) 

0.0 

0.0 

0.0 

40 

4.50  (44) 

8.37  (-11) 

0.0 

0.0 

0.0 

41 

3.00  (44) 

8.37  (-11) 

0.0 

0.0 

0.0 

42 

2.00  ( *4 ) 

8.37  (-11) 

0.  0 

0.0 

0.  0 
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4.  CALCULATIONS 


The  calculations  performed  in  generating  the  data  base 
for  the  prompt  and  delayed  dose  components  are  described  in 
this  section.  The  solution  techniques  outlined  previously 
are  further  explained  and  difficulties  encountered  are  noted 
and  discussed. 

4.1  SOURCE  ENERGY  IMPORTANCE  DISTRIBUTIONS  IN  UNIFORM 
DENSITY  AIR  OVER  GROUND 

4.1.1  Geometry  Model 

The  two  dimensional  discrete  ordinates  multigroup  code 
DOTSAI  was  used  to  calculate  the  energy  and  space  dependent 
source  importance  distributions  for  the  free- field  and  armor 
shielded  tissue  dose  detectors.  The  calculations  were  per- 
formed in  R-Z  geometry  using  the  model  depicted  in  Figure  7. 
The  point  detector  response  function  (the  adjoint  source) 
is  located  1.75  meters  above  the  ground  of  thickness  1.0 
meter.  The  model  extends  to  1500  meters  in  the  radial  direc- 
tion and  for  a nominal  distance  of  1275  meters  above  the 

- 3 3 

ground.  An  air  density  of  1.11  x 10  gm/cm  was  used. 

4.1.2  Boundary  Conditions 

The  calculations  were  performed  in  cylindrical  geom- 
etry and  the  left  boundary  is  treated  as  the  axis  of  sym- 
metry with  a reflected  boundary  condition.  The  top,  right 
and  bottom  boundaries  are  defined  with  an  albedo  boundary 
condition.  The  albedo  boundary  condition  implemented  in 
DOTSAI  is  a modified  white  boundary  condition  in  that  a 
fraction  (the  "albedo")  of  the  exiting  particles  are  returned 
with  an  isotropic  angular  distribution.  The  required  albedos 
can  be  calculated  from  one-dimensional  ANISN  calculations 
using  the  output  boundary  currents  to  calculate  the  albedo 
as  the  incoming  current  divided  by  the  exiting  current. 
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Calculated  in  this  manner,  the  change  in  energy  distribution 
of  the  incoming  particles  is  implicit  and  is  consistent  with 
the  manner  in  which  the  albedo  is  used  in  DOTSAI . For  each 
of  the  five  response  function  sources,  the  albedo  for  air 
and  ground  was  calculated.  The  differences  among  the  free- 
field  neutron,  armor  shielded  neutron  and  armor  secondary 
gamma  ray  source  albedos  were  small  as  were  the  differences 
between  the  free-field  gamma  ray  and  armor  shielded  gamma  ray 
sources.  Therefore,  only  the  free-field  neutron  and  gamma  ray 
albedos  for  air  and  ground  are  given  in  Table  12.  However, 
since  the  albedos  were  generated  for  all  given  source  spectra, 
primarily  to  investigate  any  differences,  the  appropriate 
albedo  was  used  in  the  two  dimensional  calculations. 

4.1.3  Source  Treatment 

A well  known  difficulty  in  performing  two  dimensional 
discrete  ordinates  calculations  in  a low  scattering  medium 
is  the  discrete  ray  streaming  effect.  This  effect  is  due  to 
the  preferential  transport  of  particles  along  the  discrete 
angles  of  the  angular  quadrature.  An  effective  technique  to 
help  alleviate  the  problem  is  to  analytically  calculate  the 
energy  and  space  dependent  uncollided  flux  and  generate  the 
source  particle  distribution  which  results  from  the  first 
collision.  This  technique  of  treating  the  source,  known  as 
a first  collision  source,  was  used  to  calculate  the  source 
importance  distributions. 

To  distinguish  between  the  dose  from  neutrons  and  the 
dose  from  gamma  rays  it  is  necessary  to  separate  these  por- 
tions of  the  detector  response  functions  for  use  as  adjoint 
source  terms.  Therefore,  five  calculations  were  required 
to  determine  the  source  importance  distributions  for  the  dose 
quantities  desired.  A complication  apparently  unique  to  the 
solution  of  the  adjoint  form  of  the  transport  equation  for 
the  calculations  performed  was  encountered. 
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lable  12.  Adjoint  air  and  ground  albedo  for  free-field 

neutron  and  gamma  ray  tissue  response  function 
source . 


— 

E . . (MrV) 
High 

[ N’eut  rons ) 

Neutron  Response  Source 

Gamma  Hay  Response  Source 

Air 

Ground 

Air 

Ground 

1.50  M) 

1.  568  (-1) 

1. 688  (-1) 

2.  899  (-1) 

4.734  (-1) 

1.22  Ml 

1. 771  (-1) 

2.024  (-1) 

2.934  (-1) 

5.007  (-1) 

1.00  (-1) 

1.414  (-1) 

1.930  (-11 

2.436  (-1) 

4.  934  (-1) 

8.  13  (-0) 

1.  548  (-11 

2.  060  (-11 

2 453  (-11 

5.  189  (-1) 

6.37 

1. 515  (-1) 

1.978  (-11 

2.  682  (-1) 

5.410  (-1  i 

4.96 

1.  735  (-1) 

2.  162  (-11 

3.  705  (-1) 

5.  684  (-11 

4.06 

4.283  (-1) 

3.985  (-1) 

5.904  (-11 

6.  745  (-1) 

3.01 

1. 659  (-1) 

1.773  (-1) 

4.  205  (-11 

5.  509  (-11 

2.38 

1.  034  (-1) 

4.  555  (-2) 

3.435  (-1) 

4.492  (-11 

2.30 

1.  633  (-1) 

2 446  (-11 

4.263  (-1) 

6.319  (-1) 

1.83 

2.  816  (-1) 

3.  626  (-1) 

5.285  (-11 

6.939  (-1) 

1.11  (*0) 

3.026  (-1) 

5.329  (-1) 

5.611  (-1) 

7.  560  (-1) 

5.50  (-1) 

4.466  (-1) 

6.  876  (-11 

7.  198  (-1) 

8.098  (-1) 

1.  11  (-D 

4.582  (-1) 

6.  819  (-1) 

7.983  (-1) 

7.72J  (-1) 

2.  18  (-2) 

5.  014  (-1) 

6.  626  (-1) 

8 436  (-1) 

7.369  (-1) 

3.35  (-3) 

4.957  (-1) 

6.734  (-1) 

8.  616  (-1) 

7.377  (-1) 

5.  83  (-4) 

4.  863  (-11 

6.  606  (-1) 

8.  725  (-1) 

7.242  (-1) 

1.01  (-4) 

4.407  (-1) 

6.570  (-1) 

8.  777  (-1) 

7.  183  (-1) 

2.90  (-5' 

4.207  (-1) 

6.  555  (-1 1 

8.  783  (-1) 

7.  145  (-11 

1.07  (-5) 

4.  162  (-1) 

6.553  (-1) 

8.  793  (-1) 

7.  121  (-1) 

3.06  (-6) 

3.700  (-1) 

6.557  (-1) 

8.  811  (-11 

7.  100  (-1) 

1. 12  (-6) 

3.617  (-1) 

6.567  (-1) 

8.  823  (-11 

7.090  (-1) 

4.14  (-7) 
(Gamma  Rays] 

4.  153  (-1) 

6.765  (-1) 

8.931  (-1) 

7.246  (-11 

1.00  (.1) 

0.0 

0.  0 

3.226  (-31 

1.  833  (-31 

8.  00  (+0) 

0.0 

0.  0 

3.  603  (-3) 

1. 972  (-31 

6.00 

0.0 

0.  0 

3.  660  (-3) 

1. 965  (-31 

5.00 

0.  0 

0.  0 

3.  822  (-3) 

2,  053  (-31 

4.00 

0.0 

0.0 

4.062  (-3) 

2.745  (-3) 

3.00 

0.0 

0.0 

3.  853  (-3) 

2.  143  (-3i 

2.  50 

0.0 

0.0 

3.  847  (-31 

2.  196  (-31 

2.00 

0.0 

0.0 

3.385  (-3) 

2.297  (-31 

1.50 

0.0 

0.0 

4.403  (-3) 

2.595  (-31 

1.00  (-0) 

0.0 

0.0 

6.093  (-31 

3.  163  (-3) 

7.  00  (-1) 

0.0 

0.  0 

1.313  (-2) 

5.  543  (-31 

4.50 

0.0 

0.0 

4.  818  (-2) 

2.  679  (-2) 

3.00 

0.  0 

0.0 

8.  184  (-2) 

8.  830  (-2) 

1.50 

0.0 

0.0 

9.720  (-21 

5 896  (-21 

1.00  (-1) 

0.0 

0.  0 

1.510  (-1) 

9.  769  (-21 

7.00  (-2) 

0.0 

0.0 

2.501  (-1) 

5.  877  (-2) 

4.50 

0.0 

0.0 

1.957  (-1) 

1.475  (-2) 

3.00 

0.0 

0.0 

9.  647  (-2) 

2.314  (-31 

2.00  (-2) 

0.0 

0.0 

1.073  (-2) 

1.  130  (-41 
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Multigroup  discrete  ordinates  solutions  to  the  adjoint 
transport  equation  are  obtained  by  inverting  energy  dependent 
input  data  with  respect  to  energy,  including  the  group-to- 
group  transfer  probabilities  of  the  cross  sections.  (Of 
course,  meaningful  results  are  obtained  only  with  appropriate 
source  and  boundary  terms.)  In  order  to  calculate  the  source 
importance  distribution  for  a neutron  dose  detector,  particles 
are  transported  from  a thermal  energy  upward  into  higher 
energy  groups.  The  mean  free  path  of  a thermal  neutron  in  1.11 
mg/cm  air  is  approximately  22  meters.  Since  there  are  no 
scattering  transfers  into  the  group,  the  flux  or  importance  is 
greatly  degraded  with  distance.  Even  though  the  uncollided 
component  is  calculated  analytically,  the  multiple  within 
group  scattering  and  the  upscattering  to  other  groups  is 
determined  by  iteration  within  the  code.  Lack  of  convergence 
and  negative  fluences  were  a major  problem.  Since  the  lowest 
energy  group  is  the  first  group  in  an  adjoint  problem  and  for 
this  problem  is  important  in  contributing  to  the  transport 
process,  it  appears  that  an  adequate  mesh  would  be  about  10 
meters  or  less.  This  kind  of  mesh  requirement  coupled  with 
numerical  difficulties  resulted  in  an  alternate  approach  being 
taken . 

To  minimize  problems  with  starting  with  a thermal  neutron 
source,  a scheme  was  used  to  calculate  the  importance  due  to 
both  neutrons  and  gamma  rays  and  to  gamma  rays  alone.  By  tak- 
ing the  difference  in  the  two  sets  of  results,  importance 
distributions  were  obtained  for  neutrons.  These  results  were 
then  validated  by  calculating  the  results  for  neutrons  only 
using  a fine  mesh  but  limiting  the  size  of  the  problem  to  a 
600  by  600  meter  grid.  There  was  good  agreement  between  the 
results  obtained  with  the  two  methods.  However,  due  to  nu- 
merical difficulties  and  "noise"  in  the  results  due  to  the 
interation  procedure,  results  for  low  energies  at  very  deep 
penetration  were  meaningless  and  were  set  to  zero. 


4.1.4  Angular  Quadrature 


An  Sg  symmetric  angular  quadrature  consisting  of  48 
angles  was  used  together  with  a Pg  legendre  expansion  of 
the  multigroup  scattering  cross  sections.  Results  obtained 
with  this  quadrature  indicated  severe  oscillations  as  a 
function  of  position.  These  oscillations  were  a factor  of 
two  or  more  and  therefore  resulted  in  a degradation  in 
quality  of  the  results.  Previous  calculations  also  indicated 
this  type  of  oscillation  but  detailed  results  were  not 
published . 

To  reduce  the  severity  of  this  effect,  each  problem 
was  restarted  for  one  inner  interation  per  group  using  an 
asymmetric  quadrature  consisting  of  104  angles.  (An 
quadrature  was  used  in  the  upward  direction  with  the  Sg 
quadrature  in  the  downward  direction.)  This  approach 
reduced  the  oscillations  to  a 10-307o  effect  for  some  energy 
groups  and  removed  them  altogether  for  other  groups.  It 
should  be  noted  that  the  most  severe  oscillations  occured 
for  those  energy  groups  having  a very  low  dominance  ratio, 
i.e.,  the  within  group  scattering  is  small  and  in  general 
the  mean  free  path  is  large.  Some  effects  of  the  discrete 
angles  still  exist  in  the  results  but  they  have  been  sig- 
nificantly reduced. 

Figure  8 shows  results  for  two  neutron  groups  when 
48  and  104  angles  are  used.  The  use  of  104  angles  signif- 
icantly reduces  the  oscillations  in  the  results  but  also 
appears  to  reduce  the  overall  intensity.  Figure  9 shows 
a somewhat  worst  case  for  gamma  rays.  Attempts  were 
initially  made  to  resolve  the  problem  with  an  quadrature 

(70  angles)  but  the  results  were  still  not  adequate.  As 
noted  in  Figure  9,  the  use  of  104  angles  removed  the  major 
oscillations.  All  of  the  data  base  source  energy  importance 
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GROUND  RANGE  (ml 


Figure  9.  Free-field  gamma  ray  dose  source  importance 
versus  ground  range  for  a source  1100  m 
above  ground. 


distributions  were  determined  by  performing  at  least  one 
iteration  with  104  angles. 

4.2  DELAYED  DOSE  DATA  BASE 

The  IDEA  code  and  the  source  energy  importance 

distributions  for  uniform  density  air  were  utilized  to 

obtain  the  delayed  dose  components  of  the  data  base.  A 

300.0  meter  ground  elevation  above  sea  level  was  chosen 

for  which  the  air  density  at  the  surface  is  approximately 
3 

1.19  mg/cm  . Therefore,  the  data  base  source  energy  impor- 
tance distributions  were  rho-R  scaled  from  an  air  density 
of  1.11  mg/cm  to  1.19  mg/cm  . The  IDEA  code  calculates 
the  prompt  and  delayed  radiation  from  one  or  more  weapon 
bursts.  In  the  transport  of  the  delayed  radiation,  the 
rising  and  expanding  fireball  is  considered.  The  varying 
air  mass  thickness  from  the  source  to  the  detector  due  to 
the  shock  wave  is  determined  at  each  time  step  with  rho-R 
scaling  used  to  determine  the  radiation  dose  at  the  detector. 

The  delayed  neutron  and  fission  product  gamma  ray  source 
spectra  are  given  in  Table  13.  These  two  spectra  were  folded 
with  the  importance  distributions  to  obtain  the  dose/unit 
source  particle  for  all  burst  heights  and  ground  ranges.  The 
eleven  dose  quantities  were  converted  to  4nR  values  as  a 
function  of  rho-R,  where  rho-R  is  the  slant  range  mass  thick- 
ness of  air  as  measured  from  a detector  located  1.75  meters 
above  the  ground  to  all  possible  source  points.  This  matrix 
was  then  utilized  in  a table  look-up  as  illustrated  in  Figure 
10.  The  figure  depicts  a detector  configuration  at  time  t 
seconds  after  a weapon  burst.  At  time  T=0.0  seconds,  a weapon 
burst  occurs  at  a distance  HB  above  the  ground.  As  the  blast 
phenomonlogy  is  treated  in  IDEA,  the  fireball  (the  delayed 
radiation  source)  rises  to  a point  HS  above  the  ground  at 
time  t.  The  mass  thickness  of  air,  pRS , from  the  source  to 
the  ground  is  determined  as  well  as  the  mass  thickness  of  air 
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Table  13.  Delayed  neutron  and  fission  product  gamma 
ray  source  spectra. 


Energy  Interval 
(MeV) 

Neutrons /Group 

Energy  Interv.  1 
(MeV) 

Gamma  Rays  Group 

1.83  - 1.11 

0.085 

5.0  - 4.0 

0. 0054 

1.11  - 0.55 

0.334 

4.0  - 3.0 

0. 0182 

0.55  - 0.111 

0.496 

3 0 - 2.5 

0. 0218 

0.111  - 0.0218 

0.0797 

2.5  - 2.0 

0.  0401 

0.0218  - 0.00335 

0. 0023 

2.0  - 1.5 

0. 0742 

0.00335  - 0.000583 

0.003 

1.5  - 1.0 

0.  133 

1.0  - 0.7 

0.  141 

0.7  - 0.45 

0.  159 

0.45  - 0.30 

0.  114 

0.30  - 0.15 

0.  157 

0.15  - 0.10 

0. 0429 

0 10  - 0.07 

0.  0411 

0.07  - 0.045 

0. 0532 

souncr  position 
AT  T I SECONDS 


Figure  10.  Source  and  detector  configuration  for 
calculation  of  delayed  dose  rate. 


2 

between  the  source  and  detector,  j RD  Tlu>  matrix  of  4R 
dose  versus  rho-R  is  then  interrogated  using  these  values 
to  determine  the  dose  at  the  detector.  The  dose  is  then 
multiplied  by  the  appropriate  source  intensity  and  divided 
by  the  integration  time  step  interval.  The  total  integration 
time  is  determined  to  be  the  time  at  which  the  maximum  dose 
rate  for  all  detector  responses  has  fallen  by  at  least  a 
factor  of  50.  This  procedure  was  used  for  all  of  the  dose 
quantities  for  all  of  the  burst  height-blast  yield  combinations 
of  interest. 

The  blast  yields  calculated  range  from  0.01  K.T  to  100.00 
KT  in  powers  of  10  with  additional  ones  at  30  and  500  KT.  The 
height  of  burst  for  which  these  six  yields  were  calculated  are 
13,  28,  60,  129,  278  and  476  meters.  In  addition,  each  yield 
was  calculated  at  1 meter  (surface  burst)  and  at  1000  meters. 
Thus,  a total  of  48  blast  yield-burst  height  combinations 
were  calculated.  Since  the  delayed  dose  scales  directly  with 
the  fission  yield  for  a given  blast  yield,  the  data  base  was 
generated  for  unity  fission/blast  yield  ratios.  The  resulting 
integrated  doses  for  each  of  the  eleven  detector  responses  for 
the  above  blast-burst  combinations  form  the  delayed  dose  data 
base . 
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5.  RESULTS 


The  results  of  adjoint  calculations  are  the  energy 
and  space  dependent  importance  distribution  which  may  be 
intepreted  as  the  contribution  to  the  dose  of  a neutron  or 
gamma  ray  source  with  a given  energy  placed  ar  a given  range 
and  height  above  the  ground.  This  data  may  be  used  with  any 
source  spectrum  at  any  height  to  determine  the  dose  at 
given  range.  Also,  the  data  may  be  used  to  determine  the 
relative  importance  as  a function  of  energy  of  a source 
spectrum.  Separate  data  are  required  for  each  type  of  response 
function. 

Source  energy  importance  distributions  as  a function  of 
ground  range  for  seven  source  heights  (or  heights  of  burst)  for 
the  free-field  and  armor  shielded  detectors  are  presented  in 
Tables  14  through  41.  For  each  height  of  burst  and  detector, 
the  importance  of  a source  particle  in  contributing  to  the 
prompt  neutron  and  gamma  ray  tissue  dose  is  given. 

Data  in  Tables  14  through  41  can  be  used  to  determine 
the  component  of  the  dose  due  to  different  source  energy 
spectra.  For  example,  in  Table  14  for  a ground  range  of  400 
meters  it  is  noted  that  the  importance  is  approximately  con- 
stant for  energies  greater  than  2.38  MeV.  This  means  that 
the  neutron  dose  is  not  sensitive  to  the  shape  of  the  source 
weapon  spectrum  for  neutron  energies  above  2.38  MeV.  It  is 
also  noted  that  for  neutrons  with  energies  less  than  20  keV, 
the  importance  is  a factor  of  500  less  than  it  is  for  a high 
energy  neutron.  Detailed  analysis  of  the  importances  can 
indicate  the  variations  in  dose  at  a given  range  that  result 
from  different  source  spectra.  For  example,  Table  14  would 
be  used  for  the  free  field  neutron  dose  component  for  a burst 
height  of  1 meter  and  Table  15  would  be  used  for  the  free-field 
secondary  gamma  ray  and  prompt  gamma  ray  component.  Similar 
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Table  15.  Free-field  gamma  tissue  dose  source  energy  importance 
vs.  ground  range  for  a burst  height  of  1.00  meter. 
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Armor  shielded  gamma  tissue  dose  source  energy  importance 
vs.  ground  range  for  a burst  height  of  1.00  meter. 
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Table  37.  Armor  shielded  gamma  tissue  dose  source  energy  importance 
vs.  ground  range  for  a burst  height  of  476.00  meters. 
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Table  39.  Free-field  gamma  tissue  dose  source  energy  importance 
vs.  ground  range  for  a burst  height  of  800.00  meters. 
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results  for  armor  shielded  doses,  including  the  component 
due  to  gamma  rays  produced  in  the  armor,  can  be  obtained 
from  Tables  16  and  17.  Table  17  indicates  that  low  energy, 
prompt  gamma  rays  less  than  100  keV  are  totally  shielded  by 
the  armor. 

In  order  to  present  values  of  prompt  dose  it  is 
necessary  to  fold  a source  spectra  with  the  importance 
distribution.  A detailed  example  is  given  in  Appendix  A. 

Three  source  spectra,  EM-1  fission^^,  low  yield  thermo- 
nuclear^'7^ and  Henre^^  are  given  in  Table  42  for  an  assumed 
intensity  of  each.  The  Henre  source  is  not  a weapon  source 
spectra  but  is  an  approximation  to  the  leakage  at  90  ' from  the 
thick  target  used  during  the  Henre  experiment.  These  normal- 
ized spectra  are  shown  in  Figures  11  and  12.  The  neutron 
spectra  are  given  in  Figure  11  in  terms  of  neutrons/MeV/source 
neutron  as  a function  of  energy  and  the  gamma  ray  spectra  are 
shown  in  Figure  12  in  terms  of  gamma  rays/MeV/source  gamma  ray 
as  a function  of  energy.  Folding  these  three  source  spectra 
with  the  total  source  energy  importance  distributions  for  each 
detector  results  in  total  dose  isodose  (total  prompt  dose)  con- 
tours as  shown  in  Figures  13  through  18.  These  contours  indicate 
the  dose  1.75  meters  above  the  ground  as  a function  of  height  of 
burst . 

The  contours  are  shown  normalized  to  rads/kiloton  of  blast 
and  thus  may  be  scaled  with  blast  yield.  The  isodose  contours 
for  the  fission  and  thermonuclear  sources  evidence  slight  ray 
effects  compared  to  those  for  the  Henre  source.  It  should  be 
noted  that  for  each  of  the  sources  there  are  two  burst  heights 
which  will  give  the  same  total  dose  for  ground  ranges  greater 
than  about  200.0  meters.  This  is  due  to  the  effect  of  the 
ground  in  reducing  the  dose. 
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Table  42.  Neutron  and  gamma  ray  source  spectra  for  folding 
with  source  energy  importance  distributions. 


rou 


| Neutrons] 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 

Total 

| Gamma  Ray  ] 


Upper  Energy 
(eV) 


1.50  (+7) 
1.22  (47) 
1.  00  (+7) 
8. 18  (+6) 

6.37  (+6) 
4.96  (+6) 

4.06  (46) 

3.01  (+6) 

2.38  (+6) 
2.30  (+6) 

1.83  (+6) 
1. 11  (+6) 
5 50  (+5) 

1.  11  (+5) 

2.  18  (+4) 
3.35  (+3) 

5.83  (+2) 

1.01  (+2) 
2.90  (+1) 

1.07  (+1) 
3.06  (+0) 
1.  12  (+0) 
4.  14  (-1) 


1.00  (+7) 
8.  00  (46) 

6.00  (+6) 

5.00  (+6) 

4.00  (+6) 

3.00  (+6) 

2.50  (+6) 

2.00  (+6) 

1.50  (+6) 
1.00  (+6) 

7.00  (+5) 

4.50  (+5) 

3.00  (f5) 

1.50  (+5) 

1.00  (+5) 

7.00  (+4) 

4.50  (+4) 

3.00  ( + 4) 

2.00  (44) 


Low  Yield 

1 Fission  Thermonuclear 


(Neutrons  kT) 

I I 3.15  (+21) 


0.0 
0.  0 

7.34  (+20) 
1.27  (+21) 


1.28  (+21) 
9.57  (+20) 
1.35  (+21) 


1. 67  (+21) 

2.00  (+21) 

1.34  (+21) 

2.49  (+21) 

4.89  (+21) 

5.32  (+21) 

3.83  (+21) 

5.58  (+21) 

9.32  (+20) 

8.78  (+20) 

7.70  (+21) 

5.92  (+21) 

1.68  (+22) 

1.50  (+22) 

1.  17  (+22) 

1.83  (+22) 

2.68  (+22) 

2.33  (+22) 

1.04  (+22) 

1.93  (+22) 

1.19  (+22) 

2.22  (+22) 

0.0 

5.48  (+22) 

0.0 

1.  66  (+22) 

0.0 

1.56  (+21) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

1.0  (+23) 

2.0  (+23) 

(Gamma  Rays/kT) 

4.8  (+18) 

0.0 

5.00  (+19) 

1.92  (+20) 

6.73  (+19) 

4.38  (+20) 

2.24  (+20) 

1.75  (+21) 

6.35  (+20) 

3.29  (+21) 

6.70  (+20) 

1.01  (+22) 

1.13  (+21) 

8.82  (+21) 

2.03  (+21) 

1.25  (+22) 

4.23  (+21) 

1.35  (+22) 

4.96  (+21) 

8.45  (+21) 

6.02  (+21) 

5.55  (+21) 

4.14  (+21) 

4.07  (+21) 

3.08  (+21) 

9.97  (+20) 

1.43  (+21) 

1.23  (+20) 

2.09  (+20) 

7.20  (+19) 

3.04  (+20) 

7.60  (+19) 

4.09  (+20) 

3.  10  (+19) 

4.09  (+20) 

3. 10  (+19) 

0.0 

0.0 

3.00  (+22) 

7.00  (+22) 

Hen  re 


7.1  (+23) 
8.0  (+22) 

4.0  (+22) 
2.98  (+22) 

3.02  (+22) 

3. 00  (+22) 

3.00  (+22) 
2.22  (+22) 
2.55  (+21) 
1.83  (+22) 

7.00  (+21) 

0.0 

0.0 
0.  0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 

1.0  (+24) 

1. 00  (+21) 
3.62  (+21) 
2.78  (+21) 

5.00  (+21) 

6.00  (+21) 
6.00  (+21) 

7.00  (+21) 
1.26  (+22) 
1.  64  (+22) 
1.  61  (+22) 
1.30  (+22) 
6.90  (+21) 

3.  12  (+21) 
1.75  (+20) 

1.03  (+20) 
1.  09  (+20) 

4.  43  (+19) 
4.43  (+19) 
0.0 

1.00  (+23) 
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Figure  11.  Normalized  neutron  spectra  as  a function 
of  energy. 
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Figure  13.  Prompt  free-field  tissue  dose  (rads/KT)  1.75  m 
above  ground  as  a function  of  burst  height  for 
a fission  source. 


Figure  14.  Prompt  free-field  tissue  dose  (rads/KT)  1.75  m 
above  ground  as  a function  of  burst  height  for 
a low  yield  thermonuclear  source. 
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Figure  15.  Prompt  free-field  tissue  dose  (rads/KT)  1.75m 
above  ground  as  a function  of  burst  height  for 


Figure  16.  Prompt  armor  shielded  tissue  dose  (rads/KT) 
1.75  m above  ground  as  a function  of  burst 
height  for  a fission  source. 
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The  data  in  Figures  13-18  are  similar  to  that  generally 
given  for  overpressure  results.  By  comparing  the  two  types  of 
curves  it  is  possible  to  determine  the  relationship  of  over- 
pressure and  dose  as  a function  of  the  height  of  burst.  It 
should  be  noted  that  for  large  yields  near  the  ground  the 
delayed  radiation  is  significant  and  will  result  in  the  dose 
being  larger  with  the  contours  turning  out  as  they  anproach  the 
ground.  For  lower  yields  or  higher  heights  of  burst,  the 
contours  are  not  affected  significantly. 

An  edit  code  was  written  to  produce  tables  of  dose  versus 
ground  range  for  a given  weapon  burst.  The  fission,  thermo- 
nuclear and  Henre  sources  were  edited  for  blast  yields  of  50, 
10.0,  30.0  and  100.0  KT  for  burst  heights  of  103,  129,  186  and 
278  meters,  respectively.  The  fission  to  blast  yield  ratios 
are  1.0,  0.5  and  0.1  for  the  fission,  thermonuclear  and  Henre 
source  respectively.  The  results  are  given  in  Tables  43-114 
For  each  source  and  height  of  burst,  the  prompt,  delayed  and 
prompt  plus  delayed  doses  are  given  for  the  free-field  and 
armor  shielded  detectors.  The  individual  components  of  the 
dose  are  shown  for  each  ground  range  together  with  the  neutron/ 
gamma  ray  dose  ratio. 

The  variation  of  total  dose  (prompt  plus  delayed)  as  a 
function  of  burst  height  is  shown  in  Tables  115-138  for  these 
three  sources  and  blast  yields.  In  general,  the  armor  shielded 
neutron  to  gamma  ray  dose  ratio  is  less  than  that  for  the  free 
field  case.  Although  the  gamma  rays  are  attenuated  more  by  the 
armor  than  the  neutrons,  the  shielded  gamma  ray  dose  is  en- 
hanced by  gamma  rays  produced  in  the  armor  by  the  neutrons. 

The  neutron  to  gamma  ray  dose  ratio  is  dependent  on  the  yield 
and  the  height  of  burst  and  varies  from  about  0.2  for  a 100 
KT  fission  source  to  about  15.0  for  a 100  KT  Henre  source. 

The  ratio  of  the  prompt  and  delayed  dose  for  the  fission 
and  low  yield  thermonuclear  sources  are  given  in  Tables  139  - 
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Figure  18.  Prompt  armor  shielded  tissue  dose  (rads/KT) 
1.75  m above  ground  as  a function  of  burst 
height  for  a Henre  source. 


146.  Similar  tables  are  not  presented  for  the  Henre  source 
since  the  delayed  dose  is  extremely  small  compared  to  the  prompt. 
It  is  noted  that  the  armor  preferentially  shields  the  delayed 
radiation  and  for  a 100  KT  fission  source  the  delayed  com- 
ponent becomes  equal  to  the  prompt  dose  for  near  surface  bursts. 
On  the  other  hand,  the  free-field  dose  from  delayed  radiation 
is  a factor  of  about  four  greater  than  the  prompt  radiation  for 
a 100  KT  fission  source  at  low  heights  of  burst.  For  sources 
with  less  f ission/total  yield,  the  delayed  radiation  is  not 
as  important. 
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6 COMPARISONS 


This  section  presents  comparisons  of  selected  results 
of  the  handbook  data  base  with  published  results  and/or  other 
calculations  1 methods.  In  addition,  results  from  sensitivity 
studies  performed  for  the  delayed  dose  component  will  be 
presented  and  discussed  as  well  as  a comparison  with  experi- 
mental data. 

6.1  PROMPT  AND  DELAYED  DOSE 

Dose  versus  ground  range  results  for  radiation  transport 

in  an  air  over  ground  geometry  from  a 12.2  - 15.0  MeV  neutron 

(19) 

source  and  a fission  source  spectra  were  compared  with 
values  from  the  handbook  data  base.  The  sources  are  located 
15.24  meters  above  the  ground  with  the  detector  0.5  meter 
above  the  ground.  These  two  sources  were  folded  with  the 
prompt  source  energy  importance  distribution  and  a comparison 
of  the  neutron  tissue  dose  is  shown  in  Figures  19  and  20  and 
a comparison  of  the  gamma  ray  tissue  dose  from  secondary  gamma 
rays  produced  in  the  air  and  ground  is  shown  in  Figures  21  and 
22.  The  primary  differences  between  the  previous  results  and 
those  presented  here  are  in  the  cross  sections  and  ground 
composition.  The  previous  dose  values  were  calculated  using 
Straker  cross  sections  while  the  handbook  data  base  utilized 
the  Young  evaluation. 

The  effect  of  the  different  cross  sections  is  most 
evident  in  the  secondary  gamma  ray  comparisons  with  the 
present  results  approximately  a factor  of  2 lower  than  the 
previous  results.  The  ground  composition  used  in  the  previous 
calculations  contained  8.5  weight  percent  water  compared  with 
3.0  weight  percent  for  the  handbook  results.  From  the  sensi- 
tivity analysis  of  ground  composition  and  the  effects  of  soil 
water  content,  the  predicted  dose  from  the  handbook  data  base 
would  be  expected  to  be  approximately  20?»  higher.  This  is 
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Figure  19.  Calculated  free-field  neutron  tissue  dose 
from  a 15.0  - 12.2  MeV  neutron  source 
15.24  m above  ground. 
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Figure  20.  Calculated  free-field  neutron  tissue 
dose  from  a fission  source  15.24  m 
above  ground. 


4irR?  DOSE  (cmJ  RADS/SOURCE  NEUTRON! 


DOSE  (cm1  RADS/SOURCE  NEUTRON) 


pR  (gm/cm*) 


Figure  22.  Calculated  free-field  secondary  gamma  ray 
tissue  dose  from  a fission  source  15.24  m 
above  ground. 


generally  true  for  the  neutron  dose  from  the  14  MeV  neutron 
source;  however,  the  dose  from  the  fission  source  is  approx- 
imately 40%  - 50%  higher  and  is  due  to  changes  in  gamma  ray 
production  in  the  ground.  The  combined  cross  section  and 
ground  composition  interaction  results  in  the  differences 
shown.  The  effect  of  the  difference  in  the  two  detector 
heights  is  assumed  to  be  small. 

MORSE  and  forward  DOTSAI  calculations  for  a low  yield 
thermonuclear  source  36.5  meters  above  ground  have  been  made 
as  part  of  the  analysis  of  terrain  effects  on  radiation. 

Figure  23  shows  a comparison  of  these  results^^  with  the 
handbook  data  base.  Calculational  geometry,  cross  sections 
and  material  compositions  are  the  same  and  as  expected  the 
agreement  is  very  good.  The  use  of  two  separate  calculational 
methods  provides  added  confidence  in  the  results. 

A comparison  with  results  predicted  with  the  IDEA  (75) 

code  and  the  handbook  data  were  also  made.  In  order  to  make 

a valid  comparison,  the  ground  correction  factors  from  ATR- 

(21) 

Version  4 were  incorporated.  Shown  in  Figure  24  are  the 
total  prompt  and  delayed  tissue  doses  versus  ground  range 
from  IDEA  and  the  handbook  data  base  for  a 10  KT  (507„  fission) 
typical  thermonuclear  weapon  burst  129  meters  above  ground 
The  IDEA  calculations  were  made  with  and  without  the  ground 
correction  factors  as  noted. 

The  agreement  of  the  prompt  dose  between  the  handbook 
and  IDEA  (75)  with  a ground  correction  is  satisfactory  and 
indicates  the  ground  correction  factors  for  neutrons  (which 
constitute  approximately  90%  of  the  dose)  are  adequate  The 
agreement  for  the  delayed  dose,  for  which  the  fission  product 
gamma  rays  dominate,  is  somewhat  worse.  This  indicates  that 
the  ground  affects  source  gamma  rays  differently  than  is 
presently  modeled. 
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Figure  23.  Calculated  total  free-field  tissue  dose 
from  a low  yield  thermonuclear  source 
36.5  m above  ground. 


DOSE  (RADS) 


Comparisons  have  also  been  made  with  results  predicted 

by  a simplified  tissue  dose  model  (TDM)  for  this  same  weapon 

(22) 

burst  configuration  The  total  prompt  and  delayed  free- 

field  tissue  doses  as  predicted  by  TDM  and  the  handbook  data 
base  are  compared  in  Figure  25. 

The  prompt  dose  from  the  simplified  model  is  based  on 
one-dimensional  source  energy  importance  distributions  cal- 
culated for  infinite  uniform  air.  These  importances  were 
calculated  using  the  same  cross  sections  and  air  composition 
as  were  utilized  for  the  handbook  data  base.  Correction 
factors  for  the  effect  of  the  ground  are  those  taken  from 
ATR-Version  3.  As  is  evident  from  the  comparison  shown,  the 
agreement  is  very  good. 

The  delayed  radiation  tissue  dose  model  in  TDM  is  based 
on  a graphical  procedure  to  calculate  the  fission  product 
gamma  ray  dose.  Although  a hydrodynamic  enhancement  factor 
is  utilized  to  account  for  the  disturbed  atmosphere,  no 
correction  is  made  for  the  effect  of  the  ground  on  the  delayed 
radiation.  As  is  noted  in  the  figure,  the  TDM  results  are 
approximately  507o  higher  than  those  predicted  using  the  hand- 
book data  base  and  provide  further  evidence  that  the  effect 
of  ground  on  source  gamma  rays  (based  on  older  first- last 
collision  models)  is  not  treated  properly  at  present. 

A comparison  of  the  experimental  and  calculated  fission 
product  gamma  ray  dose  rate  for  a typical  thermonuclear  source 
located  450  meters  above  ground  is  shown  in  Figure  26.  The 
detector  slant  range  is  1440  meters  from  the  burst.  The  cal- 
culated values  shown  were  obtained  using  IDEA  (75)  and  the 
methodology  utilized  to  calculate  the  handbook  delayed  data 
base.  For  times  greater  than  about  0.5  second  the  agreement 
is  fairly  good  both  in  shape  and  magnitude,  with  the  exception 
of  times  around  1.5  seconds.  A sensitivity  analysis  indicates 
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Figure  25.  Calculated  prompt  and  delayed  total  free- 
field  tissue  dose  for  a 10  KT  low  yield 
thermonuclear  burst  129  m above  ground. 
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that  the  disagreement  for  times  less  than  0.5  second  is 
partially  attributable  to  the  fission  product  gamma  ray 
source  model  utilized. 


As  noted  previously,  an  edit  code  was  written  to  pro- 
cess the  handbook  data  base.  Since  the  delayed  dose  does 
not  scale  directly  with  yield  or  height  of  burst,  inter- 
polation procedures  are  used  in  editing  portions  of  the 
handbook  delayed  dose  results.  As  a check  of  the  edit  pro- 
cedures used  to  determine  dose  values  at  intermediate  points, 
two  calculations  were  performed  using  IDEA  and  the  handbook 
methodology  for  comparision  with  interpolated  dose  values 
predicted  using  the  edit  procedures.  The  free -field  and 
armor  shielded  fission  product  gamma  ray  dose  versus  ground 
range  from  a 5 KT  burst  at  100  meters  above  the  ground  is 
shown  in  Figure  27.  A similar  comparison  for  a 200  KT  burst 
350  meters  above  the  ground  is  shown  in  Figure  28.  The 
interpolation  scheme  employed  gives  dose  values  within  about 
107o  of  the  calculated  value  and  is  therefore  adequate. 
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Figure  27.  Calculated  free-field  and  armor  shielded 
fission  product  gamma  ray  dose  from  a 
5 KT  burst  100  m above  ground. 
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6.2  DELAYED  DOSE  SENSITIVITY  ANALYSIS 


An  analysis  of  the  sensitivity  of  the  delayed  dose 
component  of  the  handbook  data  to  the  fission  product  gamma 
ray  source  model,  shock  wave  - fireball  phenomenology  model 
and  ground  activation  was  performed. 

6.2.1  Fission  Product  Gamma  Ray  Source  Model 

The  delayed  neutron  and  fission  product  gamma  ray  source 

models  incorporated  in  IDEA  are  based  on  the  assumption  that 

the  energy  and  temporal  distributions  are  separable.  A recent 

evaluation  of  the  gamma  rays  emitted  from  fission  has  been 

(23) 

performed  at  Los  Alamos'  . The  time  integrated  spectra  are 
compared  in  Table  147.  The  time  dependence  of  the  U-235 
fission  product  gamma  rays  is  shown  in  Figure  29.  This  time 
and  energy  dependent  gamma  ray  source  was  incorporated  in 
IDEA  and  several  calculations  made  to  determine  the  sensitivity 
of  the  dose  to  the  source  model.  These  calculations  did  not 
utilize  the  weapon  handbook  data  base  since  the  source  spectra 
is  not  constant  in  time  and  integration  over  the  source  energy 
importance  would  need  to  be  performed  for  each  time  period. 
Nevertheless,  conclusions  based  on  the  results  should  be 
applicable  to  the  handbook  data  base. 

The  integrated  fission  product  gamma  ray  tissue  dose 
versus  ground  range  for  a 10  KT  burst  at  129  meters  above  the 
ground  is  shown  in  Figure  30.  The  detector  is  located  1.73 
meters  above  the  ground.  Both  the  Pu-239  and  U-235  time  and 
energy  dependent  fission  product  gamma  ray  sources  were  used 
and  are  compared  with  the  time-dependent  source  model  in  IDEA. 
For  small  ground  ranges,  the  newer  source  spectra  give  a dose 
a factor  of  about  2.5  times  higher  for  both  the  Pu-239  and 
U-235  sources.  At  a ground  range  of  1200  meters,  the  dose 
from  the  U-235  source  is  507„  greater  relative  to  IDEA  while 
the  dose  from  the  PU-239  source  is  257„  greater.  The  dose  rate 
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Figure  28.  Calculated  free-field  and  armor  shielded 
fission  product  gammaray  dose  from  a 
200  KT  burst  350  m above  ground. 


Table  117  Comparison  ol  IDEA  1'ission  product,  gamma 
ray  spectrum  with  Los  Alamos  evaluation. 


Upper  Enerity 
(MeV) 

Uranium-235 

.... 

Plutonium-239 

IDEA 

LASL 

LASL 

10.  0 

0.  0 

0.0 

0.0 

8.  0 

0.  0 

0.  0 

0.0 

6.5 

0.0 

0.  0069 

0.  004 

5.0 

0.  0054 

0.  0048 

0.  003 

4.0 

0.0182 

0. 0213 

0.017 

3 0 

0.  0218 

0.  0200 

0.  018 

2.5 

0.0401 

0.  0300 

0.  027 

2.0 

0.  0424 

0.  036 

0.  032 

1.66 

0.  0644 

0.  053 

0.049 

1.33 

0.  100 

0.  089 

0.  082 

1.0 

0.  0835 

0.070 

0.  066 

0.  8 

0. 1090 

0.095 

0.  092 

0.  6 

0. 1380 

0.  140 

0.  146 

0.4 

0.  0835 

0.  084 

0.  093 

0.3 

0.  0942 

0.  143 

0.  129 

0.2 

0.  1052 

0.  141 

0.  162 

0.  1 

0.  0580 

0.033 

0.  040 

0.05 

0.  0363 

0.033 

0.  040 

ELOW  . 02 

Total  KT 

5.29(^23) 

7.  09(+23) 

■w  - — 

6.  87(f23) 
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at  a ground  range  of  500  meters  from  these  calculations  is 
shown  in  Figure  31.  Since  the  shape  of  the  dose  rate  curve 
depends  on  the  shock  wave  propagation  and  fireball  rise  rates, 
the  major  differences  are  due  to  the  changing  energy  spectra 
with  time.  For  times  greater  than  0.15  seconds  the  variations 
are  not  large  and  indicate  that  the  time-energy  separation  is 
adequate  but  possible  changes  in  intensity  may  be  required. 

6.2.2  Fireball-Shockwave  Phenomenology  Model 

The  sensitivity  ot  tne  integrated  fission  product  gamma 
ray  dose  to  the  fireball  rise  model  is  shown  in  Figure  32. 

The  three  calculations  are  for  a 10  KT  burst  at  60  meters 
above  ground  and  utilized  the  calculational  methodology  of 
the  delayed  component  of  the  handbook  data  base.  The  three 
fireball  rise  models  are  those  of  the  1974  and  1975  versions 
of  IDEA  together  with  that  of  the  RANC  code.  For  all  ground 
ranges  shown,  the  IDEA  (74)  fireball  rise  model  results  in 
the  highest  integrated  dose,  the  IDEA  (75)  model  results  in 
the  lowest  integrated  dose,  with  the  RANC  model  in  between. 

The  difference  between  the  two  IDEA  calculations  varies  from 
4.5  at  100.0  meters  ground  range  to  2.0  at  1200  meters  ground 
range . 

The  explanation  for  the  wide  variation  in  the  integrated 
dose  is  given  in  Figure  33  in  which  the  fireball  position  for 
the  IDEA  (74)  and  IDEA  (75)  fireball  rise  models  is  shown  as 
a function  of  time.  At  2 seconds  after  the  burst,  the  IDEA 
(75)  fireball  has  risen  approximately  700  meters  versus  50 
meters  for  the  IDEA  (74)  fireball.  However,  at  later  times 
the  IDEA  (75)  fireball  is  rising  slower  than  the  IDEA  (74) 
fireball.  The  complexity  of  the  interplay  between  the  ground 
and  the  rising  expanding  fireball  is  illustrated  in  Figure  34 
in  which  the  fireball  position  as  a function  of  time  for  a 
10  KT  burst  at  several  burst  heights  is  shown.  The  data  for 
the  curves  are  from  IDEA  (75)  and  are  taken  from  a portion  of 
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the  calculations  performed  for  the  delayed  dose  component 
of  the  handbook  data  base. 

The  results  shown  in  Figure  30  are  not  necessarily 
typical  since  the  modeling  of  fireball  rise  depends  on  the 
yield  and  the  height  of  burst.  In  addition,  the  delayed 
radiation  dose  depends  on  the  mass  thickness  of  air  between 
the  instantaneous  position  of  the  fireball  and  the  detector. 
Thus,  the  combined  effects  of  fireball  rise  and  shock  wave 
reflection  from  the  ground  produce  the  variation  in  results. 

The  emphasis  in  fireball  rise  modeling  has  generally  been 
placed  on  predicting  experimental  data  at  larger  times 
('v  minutes)  with  less  emphasis  on  the  early  times  seconds) 
important  in  this  problem.  Consequently,  the  delayed  radiation 
data  base  may  require  modification  if  significant  changes  are 
made  to  the  early  time  modeling  of  the  fireball  rise  rate. 

6.2.3  Effect  of  Ground  Activation 

An  investigation  of  the  dose  due  to  activated  soil  carried 
aloft  in  the  fireball  was  performed.  A 10  KT  fission  source 
for  a burst  height  of  1 meter  was  investigated.  It  was  assumed 
that  half  the  neutrons  entered  the  ground  and  produced  the 
resulting  activated  soil.  It  was  also  assumed  that  all  the  ac- 
tivated soil  was  vaporized  and  carried  aloft  in  the  fireball. 

The  IDEA  code  was  modified  to  treat  the  decay  of  the  activated 
sources  in  place  of  the  fission  product  source.  Figure  35 
shows  a comparison  of  the  free-field  tissue  dose  versus  range 
for  the  fission  product  dose  and  the  ground  activation  dose. 

It  is  noted  that  for  this  case  the  ground  activation  component 
is  insignificant. 


236 


200 


GROUND  RANGE  (m> 


COMPARISON  OF  FREE  FIELD  TISSUE  DOSE  FROM  FISSION  PRODUCT  GAMMA  RAYS 
AMO  GROUND  ACTIVATION  FOR  A 1 kT  BURST  1 m ABOVE  GROUND 

Figure  35.  Comparison  of  free-field  tissue  dose  from 
fission  product  gamma  rays  and  ground 
activation  for  a 1 KT  burst  1 m above 


7. 


SUMMARY 


This  report  documents  the  results  of  a study  performed 
to  generate  a data  base  for  predicting  the  radiation  environ- 
ment of  tactical  nuclear  weapons  for  air  bursts  over  flat 
ground.  The  radiation  environment  consists  of  the  prompt  and 
delayed  dose  for  a free-field  and  an  armor  shielded  tissue 
detector.  The  state-of-the-art  solution  techniques  and  data 
used  in  the  calculations  have  been  described  in  detail. 
Detailed  results  for  a few  typical  weapon  spectra  have  been 
presented  in  various  forms.  The  sensitivity  of  the  prompt 
and  delayed  components  of  the  data  base  to  several  important 
parameters  has  been  investigated  and  results  shown. 


Portions  of  the  total  free-field  dose  results  previously 
presented  in  Tables  115  - 138  are  sumarized  in  Figures  36  - 44. 
Shown  in  each  figure  is  the  total  (prompt  plus  delayed)  tissue 
dose  versus  ground  range  for  four  heights  of  burst  from  1 
meter  to  476  meters  above  ground  for  a given  weapon  source  and 
blast  yield.  Figures  36-39  are  for  the  EM-1  fission  source 
given  in  Table  42  for  blast  yields  of  5.0,  10.0,  30.0  and 
100.0  KT,  respectively.  The  assumed  fission/blast  yield 
ratio  is  1.0  for  the  fission  source.  Figures  40-43  show  the 
predicted  total  dose  from  the  low  yield  thermonuclear  weapon 
source  for  the  above  blast  yields.  The  fission/blast  yield 
ratio  is  assumed  to  be  0.5  for  this  weapon.  Since  the  delayed 
dose  component  for  the  Henre  source  is  insignificant,  the 
total  dose/KI  as  a function  of  burst  height  is  given  in 
Figure  44  for  this  source. 


In  general,  the  predicted  dose  in  an  a 

geometry  from  prompt  neutrons  is  in  reasonab 

previous  calculations  and  do  not  exhibit  the 

(24) 

problems  encountered  in  other  studies 

from  secondary  gamma  rays  is  about  a factor 

previously  published  data.  These  observed  d 
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Figure  36.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.75  m above  ground  as  a 
function  of  burst  height  for  a 5 KT  fission 
source . 
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Figure  38.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.75  m above  ground  as  a 
function  of  burst  height  for  a 30  KT 
fission  source. 


m 

V 

fr 


241 


DOSE  (RADS) 


GROUND  RANGE  Iml 


Figure  40.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.75  m above  ground  as  a 
function  of  burst  height  for  a 5 KT  low 
yield  thermonuclear  source. 
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Figure  41.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.75  m above  ground  as  a 
function  of  burst  height  for  a 10  KT  low 
yield  thermonuclear  source. 
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Figure  43.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.73  m above  ground  as  a 
function  of  burst  height  for  a 100  KT  low 
yield  thermonuclear  source. 
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Figure  44.  Calculated  prompt  plus  delayed  free-field 
tissue  dose  at  1.75  m above  ground  as  a 
function  of  burst  height  for  the  Henre  source 
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primarily  the  result  of  improvements  in  the  cross  section 
data  base  utilized  and  in  the  application  of  state-of-the-art 
transport  methods.  The  prediction  of  the  dose  from  prompt 
source  gamma  rays  in  an  air  over  ground  geometry  is  presented 
for  the  first  time.  The  main  component  of  the  delayed  dose 
is  that  due  to  fission  product  gamma  rays.  Uncertainties  in 
these  results  for  delayed  radiation  are  significantly  larger  than 
uncertainties  in  the  prompt  dose  due  to  several  effects.  For 
example,  the  fission  product  gamma  ray  intensity  and  spectra 
contribute  an  uncertainty  of  about  507,  with  newer  evaluations 
resulting  in  a higher  dose.  The  largest  uncertainty  in  the 
fission  product  gamma  ray  dose  is  that  associated  with  the 
fireball  rise  model.  The  calculated  dose  may  vary  from  707, 
to  a factor  of  3 depending  on  different  fireball  rise  models 
for  some  yield  and  heights  of  burst  as  a result  of  this  un- 
certainty . 

Uncertainties  in  the  predicted  prompt  and  delayed  dose 
arise  from  uncertainties  in  cross  sections,  numerical  and 
calculational  difficulties,  and  uncertainties  in  phenomenology 
models.  Cross  section  uncertainties  have  been  noted  by 
comparing  calculations  with  integral  experimental  data.  In 
some  instances  numerical  problems  arose  due  to  ray  effects 
and  iteration  convergence.  However,  in  other  cases  problems 
were  avoided  by  using  albedo  boundary  conditions.  Uncertainties 
in  phenomenology  were  investigated  for  a few  cases  to  indicate 
where  emphasis  should  be  placed  on  model  verification.  From  the 
analysis  performed  for  this  work  and  through  experience  in 
performing  similar  calculations,  an  uncertainty  table  is  sug- 
gested. It  is  realized  that  source  spectra,  height  of  burst, 
etc.  can  have  a significant  effect  on  the  uncertainty  and 
therefore  there  is  no  one  set  of  numbers  that  apply  to  all  cases. 
The  uncertainty  includes  the  effects  of  calculational  method, 
data  and  model  and  should  apply  to  the  typical  problems  involving 
yields  less  than  100  KT. 
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The  uncertainties  given  in  Table  148  are  smaller  for 
scaled  burst  heights  greater  than  200  meters  and  are  perhaps 
larger  for  lower  heights  of  burst.  Reduction  of  the  uncer- 
tainties in  the  prompt  radiation  would  require  high  accuracy 
experimental  transport  results  for  both  neutrons  and  gamma 
rays.  To  reduce  the  uncertainty  in  the  delayed  radiation, 
further  analysis  of  existing  experimental  data  and  critical 
review  of  existing  phenomenology  modeling  is  required. 


Further  observations  concerning  the  results  depend  on 
the  burst  conditions  investigated.  The  delayed  component 
of  the  dose  may  be  a significant  fraction  of  the  total  dose 
for  near  surface  bursts  with  a large  fission  yield.  In  some 
instances  the  delayed  dose  can  be  a factor  of  2 - 4 greater 
than  the  prompt.  For  the  simplified  armor  shielded  tissue 
detector  model,  the  armor  reduces  the  dose  from  source  neutrons 
and  gamma  rays  but  becomes  a source  of  secondary  gamma  rays. 

The  net  reduction  in  the  dose  due  to  the  armor  as  compared  to 
the  free-field  tissue  detector  is  not  large  and  primarily 
changes  the  neutron  to  gamma  ray  dose  ratio 
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ical  nuclear  weapons. 


been 
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Table  148.  Uncertainties 
air/ground  ge 


Dose  Component 
Prompt  Neutron 

Secondary  Gamma  Ray 

Prompt  Gamma  Ray 

Delayed  Gamma  Ray 


for  tissue  dose  in  an 
:>metry  . 

Uncertainty 
At  1000  m 
+ 25% 

+ 30% 

+ 10% 

+ 300% 

- 50% 


At  2000  m 


25 


♦ 


8.  REFERENCES 


1.  Gritzner,  M L. . "A  User's  Manual  for  the  Two-Dimensional 
Discrete  Ordinates  Code  DOTSAI",  Science  Applications , 

Inc.  . SAI - 75-  747-HU  ( 1975) 

2.  Straker,  E A.,  and  Allen,  L C.,  "A  User's  Manual  for 
the  IDEA  Code",  Science  Applications,  Inc.,  SAI - 76- 532-HU 
( 1975)  . 

3.  Straker,  E.  A , Allen,  L.  C. , and  Graham,  K.  D , "IDEA- 
An  Integrated  Dose  Environment  Analysis  Code  for  Nuclear 
Weapons  Effects",  Science  Applications,  Inc.,  To  be 
published. 

4.  Straker,  E A.,  "Investigation  of  the  Adequacy  of  Nitrogen 
Cross-Section  Sets:  Comparison  of  Neutron  and  Secondary 
Gamma  Ray  Transport  Calculations  with  Integral  Experi- 
ments", Oak  Ridge  National  Laboratory,  ORNL-TM-3768  (1972). 

5.  Straker,  E.  A.,  and  Gritzner,  M.  L. , "Neutron  and 
Secondary  Gamma  Ray  Transport  in  Infinite  Homogeneous 
Air",  Oak  Ridge  National  Laboratory,  ORNL-4464  (1969) 

6.  Young,  P.  G.  and  Foster,  D.  G. , Jr.,  "An  Evaluation  of 
the  Neutron  and  Gamma  Ray  Production  Cross  Sections  for 
Nitrogen",  Los  Alamos  Scientific  Laboratory,  LA-4725 
(1972) . 

7.  Harris,  L.,  Jr.,  Trimbell,  G.  D. , and  Young,  J.  C , 
"Integral  Experiment  to  Test  Carbon  and  Nitrogen  Cross 
Sections",  Intelcom  Rad  Technology,  DNA  2986F  (1972) 

8.  Weinert,  M. , "Omnidirektionale  Messung  des  Neutronen 
and  Gamma  Spektren  in  MeV  - Bereich  bein  Ausbreitung 
von  14  MeV  - Neutronen  in  Flussigen  Stickstoff", 
Dissertation,  Christian-Albrechts  Universitat  (1973) 

9.  Straker,  E A.,  Scott,  W.  H. , Jr.,  and  Byrn,  N.  R., 

"The  MORSE  Code  with  Combinatorial  Geometry",  Science 
Applications,  Inc.,  SAI - 72- 51 1 -LJ  (1972). 

10.  Harris,  R.  J.,  Jr.,  "Neutron  Spectrum  Measurement  in 
Liquid  Nitrogen",  Gulf  General  Atomic,  GA-9886  (1969). 

11.  Reynolds,  G.  M. , "Gamma  Ray  Transport  in  Liquid  Nitro- 
gen", Gulf  General  Atomic,  GA-9429  (1969). 


251 


12. 


Engle.  W.  W. . Jr. . "A  User's  Manual  for  ANISN,  A One- 
Dimensional  Discrete  Ordinates  Transport  Code  with 
Ansiotropic  Scattering".  Union  Carbide  Corp  , Nuclear 
Division,  K-1693  (1967). 

13.  Bartine,  D.E.,  Mynatt,  F R.,  and  Oblow,  E.  M. , 

"SWANLAKE:  A Computer  Code  Utilizing  ANISN  Radiation 
Transport  Calculations  for  Cross  Sections  Sensitivity 
Analysis",  Oak  Ridge  National  Laboratory,  ORNL-TM-3809 
( 1973)  . 

1A.  Dowler,  T.W.,  Los  Alamos  Scientific  Laboratory  (Private 
Communication  1975) . 

15.  Abbott,  L.  S.,  Claiborne,  H D. , and  Clifford,  C.  E. 
(Editors),  "Weapons  Radition  Shielding  Handbook",  DASA 
1892- A. 

16.  Huszar,  L.  and  Woolson,  W.  A.,  "Version  3 of  ATR" , 

Science  Applications,  Inc.,  SAI-75-585-LJ , (1975) 

\ 

17.  Auxier,  J.  A. , et  al,  "Nuclear  Weapons  Free  Field  Envir- 
onment Recommended  for  Initial  Radiation  Shielding 
Calculations,"  Oak  Ridge  National  Laboratory,  ORN-TM-3396 
(1972)  . 

18.  Burson , Z.,  EG&G  (Private  Communication  1970). 

19.  Straker,  E.  A.,  "Time  Dependent  Neutron  and  Secondary 
Gamma  Ray  Transport  in  an  Air-Over-Ground  Geometry", 

Oak  Ridge  National  Laboratory,  ORNL- A289-Vol . 2 (1968). 

20.  Albert,  T.  E.,  Science  Applications,  Inc.,  (Private 
Communication  1975) . 

21.  Huszar,  L.  and  Woolson,  W.  A.,  "Version  A of  ATR", 

Science  Applications,  Inc.,  To  be  published. 

22.  Albert,  T.  E.,  Graham,  K.  D. , and  Smith,  H.  T. . "Develop- 
ment of  a Tissue  Dose  Model  (TDM)  for  Weapon  Effects 
Studies",  Science  Applications,  Inc.,  Internal  Memo  (197A) 

23.  Young,  P.  G. , Los  Alamos  Scientific  Laboratory  (Private 
Communication,  1975). 

2A.  Dowler,  T.  W. , Los  Alamos  Scientific  Laboratory  (Private 
Communication,  1973). 


252 


■ \ -M  . • 


APPENDIX  A 


An  example  is  provided  Co  demonstrate  the  use  of  the 
energy  dependent  importance  function  data  such  as  that  found 
in  Tables  14  - 41.  As  noted  in  the  text,  the  importance 
function  when  multiplied  by  a source  results  in  the  dose. 

Consider  the  problem  of  calculating  the  free  field 
neutron  tissue  dose  at  a ground  range  of  800  meters  from 
a 10  KT  Henre  source  at  129  meters  above  ground.  Table  42 
gives  the  source  energy  spectrum  and  the  importance  is  taken 
from  Table  26.  Table  A1  illustrates  the  required  calculation, 
namely 


where 

D(r)  is  the  dose  at  ground  range  r, 

W is  the  total  yield  in  kilotons, 

is  the  source  intensity  in  energy  group  i, 
I^(r)  is  the  neutron  source  importance  for  energy 

group  i at  ground  range  r. 


The  results  shown  in  Table  A1  should  be  compared  with  the 
dose  given  in  Table  73.  The  slight  difference  is  attributable 
to  round-off  error. 


Table  Al.  Calculation  sheet  for  sample  problem. 


Group,  i 

Energy  Range 
(MeV) 

Source,  S. 

(neutrons  KT'group) 

Importance,  I.(r) 

(rads ^neutron  group) 

S.  *I.(r) 
i r 

(rads  KT) 

1 

15.0  12.2 

7. 10  (^23) 

9.  57  (-21) 

6.79  (+3) 

2 

12.2  10.0 

8.00  (.22) 

9.  74  (-21) 

7.  79  (+2) 

3 

10.0  8.18 

4.00  (.22) 

1.  08  (-20) 

4.32  (+2) 

4 

8.18  - 6.37 

2.98  (+22) 

1.  19  (-20) 

3.55  (+2) 

5 

6 37  - 4.96 

3.02  (.22) 

1. 44  (-20) 

4.35  (+2) 

6 

4.96  - 4.06 

3.00  (+22) 

1.  14  (-20) 

3.42  (+2) 

7 

4.06  - 3.01 

3.00  (+22) 

1.07  (-20) 

3.21  (+2) 

8 

3.01  - 2.38 

2.22  (.22) 

9.  84  (-21) 

2. 18  (+2) 

9 

2.38  - 2.30 

2.  55  (.22) 

8.98  (-21) 

2.29  (.2) 

2.30  - 1.83 

1. 83  (+22) 

6.35  (-21) 

1.  16  (+2) 

11 

1.83  - 1.11 

4.43  (-21) 

3.  10  (+1) 

12  - 23 

1.11  - 1.0  (-7) 

0.0 

-- 

0.  0 

Total 

(per  KT) 

1.0  (+24) 

1 . 00  (+4 ) 

Total  for 

10  KT 

1.0  (+25) 

1. 00  (.5) 
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